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ABSTRACT 


A theoretical discussion of the diffusion of divalent impurities in polar 
crystals of the NaCl-type is given with the aid of the Stasiw—Teltow 
association model. It is shown how accurate measurements of -the 
impurity diffusion coefficient can yield fundamental information about 
the existence and concentration of impurity—vacancy pairs (‘ complexes ’) 
in a simple and direct manner. Such measurements would therefore 
provide a valuable addition to the studies of ionic crystals already made 
of their ionic conductivity, dielectric loss, ete. The main restriction 
of the present theory is that the impurity—vacancy complexes should 
have a mean lifetime which is long compared with the mean jump times 
which govern their changes of position and orientation. 


§ 1. INTRODUCTION 


A KNOWLEDGE of the properties of defects in ionic solids is essential to 
the understanding of a wide range of physical and chemical phenomena 
(e.g. colour centres and oxidation of metals and alloys). The absence of 
complications such as electronic conductivity in crystals like NaCl 
and AgBr has permitted the extensive study of their lattice defects 
(ef. Seitz 1951, 1954), and the possibility of forming substitutional solutions 
with for example CaCl, has allowed the study of systems containing 
constant controlled numbers of vacancies present in much greater concen- 
trations than the intrinsic or thermal defects. This useful fact results 
from the condition of electroneutrality, which requires that every Ca?+ 
ion dissolving substitutionally in the NaCl lattice, be accompanied 
somewhere by a Nat vacancy. The solution of, say, AICI, in a similar 
way would introduce two Nat vacancies for every Al?* ion. 
Unfortunately, however, these solid solutions are not as simple as 
appears at first sight, owing to the long range Coulomb interactions 
existing between the unbalanced charges on the impurities and the 
vacancies. This system of impurities and vacancies dispersed throughout 
an otherwise perfect crystal forms a sort of electrolyte solution rather 
than an ideal mixture. For the statistical description of these systems 
one has available on the one hand the theories of * strong ’, i.e. totally 


* Communicated by the Author. 
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dissociated, electrolytes (e.g. the Debye and Hiickel (1923) theory ; 
general references are Falkenhagen (1934) and Harned and Owen (1950)), 
and on the other the ‘weak’ electrolyte theories which describe the 
state of the system in terms of the numbers of opposite ‘ charges ’ (i.e. 
impurity and vacancy) which are associated into pairs. In the case of 
divalent impurity ions in 1:1 salts these pairs are electrically neutral 
and do not contribute to the ionic conductivity*. 

The simplest form of the weak electrolyte theory, in its application 
to the present problem, was given by Stasiw and Teltow (1947); it 
neglects the long range interactions between the unassociated impurities 
and vacancies and in the limit of high temperatures goes over into the 
theory of an ideal solid solution. In effect, it treats the equilibrium 
existing between the associated pairs and the unassociated vacancies 
and impurities as an equilibrium for a reversible chemical reaction ; the 
associated pairs, or complexes as they are often called, are treated like 
dissociating diatomic molecules. Detailed applications of this approach 
have been made among others by Teltow (1949, for AgBr+CdBr,), 
Etzel and Maurer (1950, for NaCl+CdCl,), Stasiw (1952, for silver 
halides+-Ag,S) and by Haven (1953, for NaCl+ CaCl, and Mn(Cl,). 
An extension of the theory so as to include the effects of the long range 
interactions has been given and applied to the NaCl+CdCl, system 
(Lidiard 1954). 

In general one would expect the association theory to give a better 
approximation to these crystal systems than the ‘strong’, i.e. totally 
dissociated, electrolyte theories. Thus it is known that salts such as 
KCl which are ‘ strong’ electrolytes when dissolved in water (dielectric 
constant ~ 80), must be regarded as only ‘ weak’ electrolytes when 
dissolved in solvents of much lower dielectric constant. We would 
therefore expect to find a similar situation existing within the solids 
themselves since their dielectric constants are not large. 

Granted that appreciable association does occur between impurities 
and vacancies we must next consider how to obtain the degree of associa- 
tion from experimental data. This information can be obtained indirectly 
from studies of ionic conductivity (d.c.) as a function of temperature 
and impurity content. More immediate information can be obtained 
from alternating current studies since the dipolar nature of the complexes 
gives rise to a characteristic Debye-type loss proportional to their 
concentration (Haven 1953, Lidiard 1955). A third method which is 
discussed in the present paper is the study of the diffusion of the impurities 
themselves. For the case where the mean lifetime of the complexes 
is much greater than the mean times characterizing the re-orientation 


* This would not be true in the case of a pair formed from a trivalent cation, 
e.g. Fe®+ and a Nat vacancy since here the net charge is positive. Thus 
Kelting and Witt (1949) in their experiments on KCl found that Ca2+ additions 
were not removed electrolytically although an appreciable fraction of the back- 
ground impurity in the “pure”? KCl crystals was removed in this way (pre- 
sumably polyvalent ions). : 
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jumps, we show that measurements of impurity diffusion can provide 
direct information about the degree of association as function of impurity 
concentration and temperature. Indeed it would seem strange if this 
were not so, since the impurity ions are only in a position to move when 
there is a vacancy adjacent to them, i.e. when they are in an associated 
condition. However, when the lifetime of the complexes is relatively 
short the results are in general more complicated and appear difficult to: 
apply at the present time. i 

A previous theoretical treatment of the subject of diffusion of impurities 
in simple polar crystals by Schéne, Stasiw and Teltow (1951), we believe 
to be incorrect. Our principal objections are to their use of a diffusion 
coefficient for wnassociated impurity ions and to their neglect of the 
dynamical association—dissociation equilibrium between the complexes 
and the unassociated vacancies and impurities. As noted above it is 
obvious that at any instant the only impurity ions free to move are those 
which are associated with vacancies at that instant. Therefore provided 
we take adequate account of the dissociation—association ‘ reaction ’ and 
provided we correctly describe the motion of the complexes, we have the 
basis of a correct calculation. 

In the present paper we derive our principal results in a plausible 
elementary manner (§2) and discuss their implications (§§ 3,4). The 
dependence of the diffusion coefficient on temperature and concentration 
is shown both analytically and graphically. In a subsequent paper we 
shall supply the full rigorous kinetic derivation of the general equations. 


§ 2. Impurity DirFrusion CoEFFICIENT, D(c) 


In the following calculation we have a definite steady state situation 
in mind. We suppose that we have a mixed crystal, e.g. NaCl containing 
a mole fraction, c, of CdCl,, set between a source and a sink of the impurity 
and across which there is a small steady gradient of impurity concentra- 
tion, dc/dxz. We further imagine that our crystal is in the shape of a right 
cylinder so that the problem is one-dimensional. Corresponding to the 
concentration gradient there will be small steady current of impurity 
ions J (number of ions per cm? per sec) given by Fick’s law 


J=—D(o) = (No) ie Tusa greg oer) 


where NV, is the number of cation sites per cm®. The diffusion constant 
D(c) is a function of the crystal composition. By using the Stasiw and 
Teltow (1947) association theory we shall calculate D(c) as a function of ¢ 
for the case of divalent additions. For definiteness we shall restrict 
ourselves to crystals of the NaCl type where the cation sublattice is 
face centred cubic. 

When the molar impurity concentration is ¢ let the number of impurity 
ions per cm? be V,;; then V,=N,c. Of these let be associated with 
vacancies. Since we shall assume sufficient impurity for the thermal 
defects to be negligible it follows that there are N;—n unassociated 


3K2 
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impurity ions and V,—n unassociated vacancies. If local equilibrium 
of the association—dissociation reaction is maintained then (Lidiard 1954) 

al? exp (AG/RT),~22 ly as wanna) 
where the factor 12 comes from the orientational entropy of the 
complexes. 4G is the Gibbs free energy of association and k is Boltz- 
mann’s constant. We shall return to this equation later. 

When we now consider the impurity diffusion we see at once that on 
the vacancy diffusion model the unassociated impurity ions are unable 
to move. Also the unassociated vacancies are prevented from diffusing 
to a region of lower impurity content by the Nernst diffusion potential 
which arises from the immobility of the oppositely charged impurity 
ions (see e.g. Denbigh 1951, especially §5.3). Hence in calculating the 
impurity diffusion we need to consider only the independent, random 
motions of the complexes. 

For simplicity let us study diffusion along a [100] direction in a 
cylindrical crystal of cross section 1 cem?. (NaCl-type crystals, being 
cubic, display isotropic diffusion properties so our choice of the [100] 
direction will not restrict our results.) A single lattice plane perpen- 
dicular to the diffusion direction may then be regarded as occupying a 
volume a cm® where a is the smallest anion—cation separation. Hence 
this plane will contain an associated impurity ions, i.e. an complexes 
with the impurity end in the plane x=0. One-third of these will have the 
vacancy located in the forward neighbouring plane (w=a). At any 
instant only complexes with this orientation are in a position to contribute 
to the forward diffusion of impurities. Likewise the only complexes 
able to contribute to diffusion backwards into the plane z=0 are those 
with impurity at «=a and vacancy at «=0. 

Provided that the continued diffusion of the complexes is not hindered 
by a very slow motion of the attached vacancy around the impurity 
ion (see below) the diffusion coefficient can then be calculated very simply 
in the following manner. Let w, be the probability per unit time that 
an impurity ion will exchange places with its associated vacancy. Then 
the number of impurity ions crossing a mathematical plane set between 
the planes of atoms x=0 and x=a is a n(0)w,/3 per cm? per sec in the 
forward direction and an(a)w,/3 per cm? per sec in the backward direction. 
The net number crossing is therefore 


ae dn 
J=-Su(Z), 


where we have replaced n(a)—n(0) by a(dn/dx) evaluated at x=0. 
Introducing the degree of association p = n/N, and the fractional molar 
concentration ¢ we find that a comparison of (2.3) with (2.1) yields 


d 
D(c)=34°W2 = (pe), oh pols Tas wed rc lar nae aera) 


where p is determined as a function of ¢ by (2.2). 
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Before proceeding to discuss the dependence of D on ¢ we may 
generalize the constant factor in (2.4). Our elementary discussion has 
overlooked one important feature of the process of diffusion via 
complexes ; namely that in order to continue its forward diffusion the 
complex must go through a cycle of orientation changes. Each impurity— 
vacancy exchange must be followed by two successive forward jumps of 
the vacancy around the impurity as centre—as it were to reset the 
complex for a fresh forward jump of the impurity. Hence the probability 
occurring in (2.4) is the probability of impurity—vacancy exchange only 
if the speed of the diffusion is not limited by a slow motion of the vacancy 
around the impurity. The appropriate factor in intermediate cases can 
only be found rigorously by a full kinetic discussion ; such a calculation 
shows, however, that it is also given correctly by the following plausible 
argument. 


Fig. 1 


X = IMPURITY ION 
© = NEAREST NEIGHBOUR 


The twelve nearest cation neighbours to an impurity ion (X) in a NaCl-type 
lattice. If the impurity is part of a complex, one of these twelve sites 
will be vacant. When the diffusion direction defines the v-axis as shown, 
then a, 6 and c represent the planes x=a, x=0 and x=—a. The 
impurity ion lies at the centre of «=0. 


We notice that if w, be the probability per unit time that an associated 
vacancy will jump from one attached position to another particular 
attached position, then 2w, is the probability per unit time that an 
associated vacancy will make one of the necessary forward jumps (say 
c to 6 in fig. 1) since there are two distinct ways in which it can do so. 
Likewise 2w, is the probability for executing the second forward jump 
necessary for completion of the cycle. 
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Now in general it is possible to associate with such probabilities a mean 
time which must elapse before the event occurs. Furthermore, the time 
required is the reciprocal of the corresponding probability. Thus the 
mean time required for the impurity—vacancy exchange is w,!. Now it 
is reasonable to suppose that the mean time required for the performance 
of several steps is the sum of the mean times for the individual steps. 
Hence the mean time associated with the cycle of jumps necessary for 
the continued diffusion of the complexes will be Wo1+ (2w,) + (20,)7}. 
The corresponding mean probability is the inverse of this or ww/(w,+-W9). 
It is this mean probability which must be inserted in place of w. in 
(2.4), giving 
a WW, d ote 


We observe that when w,> ws», (2.5) goes over into (2.4) as it must do. 


INC l= 


§ 3. DEPENDENCE OF D(c) ON c ACCORDING TO THE STASIW—TELTOW 
ASSOCIATION THEORY 


Equation (2.2) can be rewritten as 


a =12¢ exp (AG/RT) = 12c exp (— AS/k) . exp (4H/RT), 
so that the dependence of D(c) upon ¢ implicit in (2.5) can be expressed by 
D(c)=Dy{1—[1+ 48c exp (4G/RT)]-¥/7}, Series Ae td . 
where Di =G7 WW fs Ws)-- eed 


D(c) depends on ¢ as shown graphically in fig. 2. We notice two 
limiting cases : 


(i) cexp (AG/RT) <1, i.e. little association, then 
D(c)=24D,cexp(4G/kRT) . .. . . (3.8) 

i.e. D(c) is proportional to c. 

(ii) c exp (AG/RT')> 1, i.e. almost total association, then 

D(c) + Dg, 

i.e. D(c) tends to a saturation value. 

Naturally the practical attainment of saturation is easier at low 
temperatures, where exp (4G/R7') is large, than at high temperatures. 

The matching of experimental isothermals to eqn. (3.1) is a matter of 
inferring the unknowns D) and 4G. This may be done by inserting 
into (3.1) corresponding values of D(c) and ¢ at two points and solving 
the resulting equations for Dy and exp (4G/RT). Alternatively if the 
experimental values at low concentrations are sufficiently accurate, in 
place of one of these equations we can use the simpler relation, 


. aD 
24D, exp (4G/RT)= Lim = (cf. eqn. (3.3)). 
c>0 
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It is clear that given diffusion isotherms at a number of temperatures 
we can infer Dy and AG as functions of temperature. Furthermore, since 
dielectric loss measurements can provide values for W,+W, and aise 
Dy immediately enables w,w,/(w,+-w,) to be found, both w, and wy, can 
be found separately. Simultaneous diffusion and alternating current 
studies can therefore in principle, yield w,, w, and AG, ie. all the relevant 
parameters for tightly bound complexes. 


1O0O0c. 


The divalent impurity diffusion coefficient D(c) as a function of the impurity 
concentration (mole fraction) for three different temperatures. The 
numbers on the curves represent absolute temperature measured in 
units of 4G/k. When 4S=0 and 4H is a constant these numbers are 
proportional to 7’. 


§ 4. AGREEMENT WITH EXPERIMENTAL RESULTS ON 
Impurity DIFFUSION 


The only experimental results known to us on the kind of system we 
have considered are those obtained by Schéne, Stasiw and Teltow (1951) 
on AgBr containing CdBr, and PbBr,. These measurements are 
unfortunately not very exact and are not susceptible to a detailed analysis. 
As may be seen from tables 1 and 2 the diffusion coefficients which have 
been obtained are averages over ranges of impurity content. Nevertheless, 
it is clear that the general trend of these values corresponds with the 
predictions of our theory, particularly the tendency to saturation at high 
impurity concentrations. If we interpret the highest D value at each 
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temperature as D, then we find that the five D, values for Cd** can be 
roughly approximated by the following function of temperature 

Dy ~ 0:21 exp (—8300/7’) cm? sec™?. ig ie (4.1) 
If we suppose that the barrier height for the jump involving the impurity 
ion is great owing to its double charge then w,<w, and (4.1) and (3.2) give 


Wy ~ 7:7 X 10" exp (— 8300/7) sec 


for Cd2+ in AgBr. However, this result is clearly very uncertain. 


Table 1. Diffusion coefficient of Cd?+ ions in AgBr in units of 
10-7 cm? sec! (after Schone, Stasiw and Teltow 1951) 


tataeeal Malo 400°C 350°e | 300% | 250° | 200% 
be ré) 


0 - 0-1 
0-691 
2 -4 


Table 2. Diffusion coefficient of Pb?* ions in AgBr in units of 
10-7 em? sec~! (after Schéne, Stasiw and Teltow 1951) 


Concentration 
Interval Mol. % 


§ 5. CONCLUSION 


We have shown how measurements of the diffusion of divalent impurities 
in alkali halides, silver halides, etc. can yield direct information about 
the existence and concentration of the impurity—vacancy complexes 
employed in the Stasiw—Teltow association model. The dependence of 
the diffusion coefficient on impurity concentration and on the funda- 
mental parameters of the complexes has been calculated for the case 
where the complexes are tightly bound, i.e. for the case where their mean 
lifetime is long compared with the mean jump times involved. This 
restriction will be removed in a subsequent paper, in which we shall 
give a full kinetic discussion. However, the results for the case where 
the complexes are only loosely bound are complicated and do not appear 
convenient for the quantitative analysis of experimental data. 
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Notes added in proof 


1. Dr. Teltow has pointed out that if the terms involving the diffusion 
coefficient of the unassociated impurity ions are removed from the Schéne, 
Stasiw and Teltow theory, the resulting equations are largely equivalent 
to the equations of the present paper. In particular their eqn. (8) becomes 
identical with eqn. (2.4) above (J. Teltow, communication to be published 
in Phil. Maq.). 

2. The completed investigation of a system of complexes of arbitrary 
lifetime shows that provided w,<w, eqn. (2.4) remains true. In this 
case the diffusion isothermals are exactly as shown in fig. 2. 
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[Received March 24, 1955] 


SUMMARY 
The cross section for the reaction !°B(p«)’Be has been measured 
absolutely for proton energies of about 200 kev using magnetic analysis. 
Values of 0-92--0-13 mb, 0:66-0:09 mb and 0-29+0-04 mb were found 
for proton energies of 205 kev, 200 kev and 180 kev respectively. Cross 


sections relative to these values were obtained for proton energies down 
to 70 kev. 


$1. INTRODUCTION 

THE absolute yield of the reaction 

B+ 1H->7Be+4He+1:147 Mev . .. . . (1) 
has been measured for proton energies between 230 kev and 1600 kev 
(Burcham and Freeman 1950, Brown, Snyder, Fowler and Lauritsen 
1951). The present paper reports an extension of these observations 
down to a proton energy of 70 kev. The interest of work such as this is 
not only that evidence may be obtained for resonance levels in the 
intermediate nucleus 14C, but also that accurate measurements of the 
cross section for (px) reactions in the 100 kev region may be useful in 
astrophysical theories (Salpeter 1953). 

For proton energies above 300 kev the alpha particles emitted at 
90° to the incident beam in the !B(pa) reaction have a range less than 
that of the scattered protons and a magnetic analyser was used in the 
earlier work to discriminate between the two groups of particles. The same 
method was adopted in the present experiment for proton energies down 
to 120 kev ; below this energy the yield was found to be too small for 
convenience in using the magnetic analyser. Fortunately the alpha 
particle range for protons of 120 kev energy is sufficiently greater than 
that of the scattered protons to allow discrimination by absorbers, and 
a simple target assembly with a detector subtending a large solid angle 
could be used. 

§ 2. EXPERIMENTAL APPARATUS 

Protons extracted from a radio-frequency ion source (Thonemann, 
Moffat, Roaf and Sanders 1948) were accelerated in a vertical two-stage 
vacuum tube. The high-tension supply was a four-stage cascade 
generator using selenium rectifiers and with an output ripple, after one 


* Present address : Windscale Works, Sellafield, Cumberland. 
+ Communicated by W. E. Burcham. 
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stage of resistance capacity filtering, of 0-3°/ under the conditions of the 
experiment. The accelerating voltage was measured with an accuracy 
of 1% by a 1300 megohm resistor chain and a microammeter, The 
protons from the accelerating tube passed through a deflecting magnet 
on to the 1B target. A stop limited the useful target size to 6x 6mm 
and a negatively charged suppressor ring immediately beneath the stop 
prevented secondary electrons reaching the target from exposed edges. 
Target currents were always measured with the target at a sufficiently 
large positive potential to retain secondary electrons. An integrator 
circuit showed when the target had received a pre-determined charge. 

The °B targets used during the experiments were prepared at the 
A.E.R.E. Harwell by electromagnetic deposition of the isotope on copper. 
Two thin targets of 90 ugm em? and 120 u.gm em? (accurate to 10%) 
and one thick target (400 .gm cm~2) were used for the cross section 
‘determinations. The thickness of this last target was determined by 
weighing the backing foil before and after deposition. 


Fig. 1 


THREADED RING 


ANALYSER 
MA CNET 


renee SPECTROMETER _MAGNET- 
callingaToR ,7— aS 


LTS) 
ea) 
sat 


TY —SSupvol INSULATION 


| ~3e0v ror + 
‘SUPPRESSOR GLASS METAL SEAL 


Magnetic analyser, target and detector assembly. 


The analyser used for discriminating between alpha particles and 
protons was a length of 3 cm x 1 cm wave guide bent into a 60° arc and 
placed between the poles of an electromagnet (fig. 1). The maximum 
field available was 6000 gauss and the mean radius of the analyser vacuum 
system was 21 cm. Straight sections of waveguide, each 11 cm long, 
extended beyond the analyser field and joined the analyser at one end 
to the target assembly and at the other end to a zinc sulphide screen and 
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photomultiplier (EMI type 6260) for particle detection. Preliminary 
experiments with polonium alpha particles of differing residual ranges 
indicated that the screens, which were about 1-5 mgm cm~? thick, could 
be relied upon to count more than 95% of the incident alpha particles. 


§ 3. DETERMINATION OF THE ABSOLUTE CROSS SECTION 


To determine the cross section for the !°B(px) reaction it was necessary 
to measure : 

(i) The total number, N,, of alpha particles emitted into the analyser 
for a given number, J, of incident protons. 

(ii) The number of disintegrable nuclei per cm? in the target. This 
is nt where n is the number of disintegrable nuclei per cm? and ¢ is the 
thickness of target from which alpha particles were received into the 
analyser. 

(iii) The solid angle 2, of the analyser in the centre of mass system 
of coordinates for the °B(pa) reaction. 

The total cross section is then 


NI 
: errenee amen as 1) 


assuming an isotropic distribution of alpha particles. 

The determination of NV, and the calculation of ¢ depend on whether 
the energy loss in the target for the alpha particles is greater or less than 
the energy interval accepted by the analyser. Snyder, Rubin, Fowler 
and Lauritsen (1950) show that for thin targets, 


o=—4r 


Ws n di 3 
oN, Qant. earn . . : . . . : ( ) 
where ¢ is the actual thickness of the target, and for thick targets 
0H, cos 6, 
a 4a Neaxk “4! 0k, 1 oR Bast co beeen 4) 
Ng Sass 2H, cos 0, 


In these formulae 
(i) V, is the number of alpha-counts observed for N,, protons when 
the analyser magnet current setting is 7 (thin target). 

(ii) R is the resolution of the analyser, equal to 2H/4H where ABE is 
the range of energies of particles accepted at the energy setting ZF. 

(iii) Nmax is the maximum number of alpha particles observed for Ne 
protons over the flat portion of the alpha particle profile from a thick 
target. 

(iv) #, is the energy of outgoing alpha particles. 

(v) 6, and @, are the angles made by the incident protons (energy E,) 
and observed alpha particles with the normal to the target surface. 

(vi) «, and e, are the stopping cross sections of B for protons of 
200 kev and alpha particles of 840 kev, taken to be 440 kev mgm~! cm? 
and 1820 kev mgm! cm® respectively (Allison and Warshaw 1953). 
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The number of alpha particles observed was in each case multiplied by 
a factor of 1-64 (Dissanaike 1953) to take account of the particles emerging 
from the target in the singly charged state. Figure 2 shows a typical 
alpha particle counting-rate profile obtained from the thick target. 
For the thin targets the integral in eqn. (3) was evaluated graphically 
from curves of alpha particle counting rate against analyser magnet 
current. 


Fig. 2 


a particle yield per proton per 47 steradians x 10° 


1-0 1.5 2-0 2-5 3-0 3.5 
Analyser magnet current in amperes 


Alpha particle profile obtained from thick 1B target (400ugm cm~-*) for 
E,=200 kev. 


In the present experiment 6,=45°, and #,=19-6°. The factor 2,/R 
was found by observations on protons scattered elastically from thick 
targets of aluminium, carbon and beryllium. The solid angle 2, in the 
centre of mass system depends on the target material in the scattering 
experiments, and was transformed to laboratory coordinates ({2;) in 
each case. The ratio Q,/R is an instrumental factor which should be 
independent of the reaction being studied. 


§ 4. RESULTS 


The values of 2,/R found from the scattering experiments are shown in 
table 1. 

Using the mean of these values and transforming 2, back to the centre 
of mass system for the !°B(p«) reaction the cross ections calculated for the 
10B(po) reaction for various targets are shown in table 2. 
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The errors are based on estimated accuracies of determination of the 
quantities in eqns. (3) and (4). From the values listed we obtain 


o[°B(pa) at 180 kev]=0-29-+0-04 
o[!°B(px) at 200 kev]=0-66-+-0-09 +} millibarns*. 
o[B(px) at 205 kev]=0-92+0-13 } 


Table 1 


Proton energy Q,/R 


Scatterer Son sy ne (laboratory system) |steradians 
ogc kev (<x 10-5) 
Be 125-2° 205 1-86 
C 124-0° - 205 1-88 
Al 121-8° 1-82 


Table 2 
‘ Proton energy Reaction cross section 
Hose ERO SEES (laboratory system) (102 cm?) 
(ugm cm~*) (kev) 
90 177 2-9 +0-4 
90 _ 199 6-6 +0-9 
90 205 7441-1 
90 205 9-7 +1-4 
120 205 9-0+1:3 
400 205 10-9+1-4 


Relative values of the cross section at lower energies were obtained. 
using a thick target and the analyser for proton energies down to 120 key. 
For energies below 120 kev a simple target chamber employing range 
discrimination, and with a much larger solid angle than the analyser 
was used. The alpha particles were detected by a zinc sulphide screen 
and photomultiplier. These observations were fitted to the absolute 
determinations at 180 and 205 kev. 

The excitation curve so found is shown in fig. 3, which also includes 
the low energy points obtained by Burcham and Freeman (1950). A 
plot of log (c#) against H-!/? (fig. 4) showed no anomalies down to 
H,=70 kev and the slope of this line was within 2° of that calculated 


from the formula 
tie A 27 Ze2 
foals AV eae 


p p 
in which v, is the velocity in the centre of mass system of the proton 
incident on B. The variation of the barrier penetrability for the 


ee eaeen ie eae Te 
* These values differ slightly from those communicated to Ajzenberg and. 
Lauritsen (Rev. Mod. Phys., %7, 77, 1955). 


Fig. 3 


Cross section cm2 


100 200 300 400 500 


Proton energy 7ev in Laboratory system 
Excitation curve for the reaction 1B(pa«)7Be. 
© from Burcham and Freeman (1950); -+present work, absolute points ; 
; e present work, relative measurements. 


Fig. 4 


Cross section x proton energy cm? x kev 


2052-40 760." 80— 100-120 


(Proton energy)-? kev-! x 10° 
Gamow plot for the reaction B(pa)7Be. 


Afrom Brown e¢ al. (1951) ; © from Burcham and Freeman (1950) 
+ e@ present work. 
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emergent alpha particles was less than 39% over the proton energy range 
investigated and has been neglected. No effects due to resonance levels 
in the 4C nucleus in the energy range investigated were found. 
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ABSTRACT 


The problem of determining the time of relaxation due to scattering 
of electrons and holes by charged impurity donors and acceptors in 
semiconductors is reconsidered in an effort to derive formulae more 
exact and general than those due to Conwell, Weisskopf, and other 
authors. 


$1. INTRODUCTION 


THE author and Dr. S. K. Roy (1955) have recently computed the chief 
integrals which occur in the theory of the electrical, thermal, thermo- 
electric, optical and magnetic-transport properties of isotropic or 
polycrystalline elemental semiconductors, on the assumption that the 
most important scattering processes are those due to (1) thermal 
agitation with scattering proportional to H1/?, where H=electron energy, 
.and (2) the presence of charged impurity ions with scattering roughly 
proportional to H~-/2, In connection with these tabulations, it seems 
desirable to re-examine the theory of the latter scattering process in an 
effort to establish formulae more accurate than those discussed in the 
literature (e.g. Conwell and Weisskopf 1946, 1950, Johnson and Lark- 
Horovitz 1947, Shockley 1950, pp. 258-264), since these largely rely on 
rather arbitrary cut-off approximations. 

Our basic approach is essentially the same as that employed by Mott 
(1936) in a discussion of the resistance of dilute solid solutions of one 
metal in another. Substantial alterations are, however, required in 
Mott’s discussion in order to make it applicable to semiconductors, and 
we shall therefore re-derive the theory, taking due account of 


(a) the dielectric constant, which must be retained at all stages, 
(b) the simultaneous presence of holes and electrons, and 
(c) the degree of degeneracy of the electron and hole systems. 


§2. THE PoTENTIAL ENERGY 
Screening of Ionized Donors by Electrons 


For the sake of simplicity in exposition, we shall first consider an 
ordinary n-type impurity semiconductor in which one excess electron 


* Communicated by the Author. 
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is provided by each donor, which therefore becomes singly positively 
charged. The first approximation to the potential energy V of an 
incident electron distant 7 from this positive charge, V=e?/«r where « is 
the dielectric constant of the bulk material, would lead to an infinite 
scattering cross section for incident electrons passing the charged impurity 
ion, on account of the high probability of small-angle scattering of 
electrons travelling past at large distances, i.e. of those electrons with 
large impact parameters. In the well-known treatment by Conwell and. 
Weisskopf (1946, 1950) this difficulty is surmounted by assuming that 
V=e2/«r for electrons approaching with impact parameters smaller than 
4N 4/3, where Ng, is the number of impurity donor ions per unit 
volume, while V=0 for larger impact parameters. Put another way, 
large ‘ misses’ are simply ignored in their treatment. In order to set 
up a more accurate theory, it is necessary to find a more comprehensive 
form for V(r), particularly when r is large, taking explicit account of 
the screening of the positively charged donor ion by the surrounding 
electron cloud. 

Let n represent the average number of electrons per unit volume in 
the conduction band (excluding of course electrons actually bound to 
donors), and ¢ the electrostatic potential. Then, since in the presence 
of the self-consistent field due to ¢ the customary space-independent 
Fermi energy ¢,, must be replaced by ¢,,+-ed(r), e being taken as positive, 
the number of electrons per unit volume near a positively charged donor 
ion will be 

8974/2 m,,8/2 pe E)2dE 
n(rj= | weal Jf ES SOPOT) GRE oe 
i J. exp[{H—L,—ep(r) MTT 1 
where m,, is the effective electronic mass, assuming the energy band to 


be of standard form (energy proportional to square of wave-vector). 
Adopting the notation 


1p 2* da 
Eines ce | eee ee 
for the Fermi—Dirac integrals,* this may be written 
4arm,,3/2( 27) 1/2(kT')3/2 ed(r 
A A ep {nt Sa where ,=¢,/kT. . (3) 


Since we are primarily interested only in the form ¢(r) takes for large r, 


* This notation differs from McDougall and Stoner’s (1938) by the presence 
of the additional factor 1/(k!), and possesses the following advantages : 


(a) In the classical limit 7 <0, ¥;,(y) is independent of the order k. 


(b) Interpolation is possible between different orders k as well as between 
different arguments 7. 


(c) ae eee between the function and its derivative is simpler (see (4) 
above). 


The required integrals ¥,,, and ¥_,/. are tabulated in this form in Dingle and. 
Roy (1955). 
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where ¢ is necessarily very small, this may be expanded, making use of 
the aes relation 

Fly (n)=F y_1(7). = Se Series (C3 
Thus, to sufficient accuracy, the space-dependence of the screening 
charge is 


Arre?m,,3/2(QrkT)1/2 
p(7)=—ebn(r)=— — F _yj(n,)P(r). . « (5) 
With this form for p(r), the solution of Poisson’s equation 
V?2b=—4rp/k, pat SP ccsiel tition we. 2a (6) 


which is spherically symmetric, vanishes for roo, and tends to the 
known value 4-e/«r for r+0, where the screening cannot be effective, is 


€ 
o(f)= > eXP (—r/B), eee ee ta CT) 
the ‘ screening distance’ R being given by 
167e?m,,3/2(27kT')1/2 
UP So earns Xi gmeremmmmce rTP, 
= 9-923 x 10-8cl2P-1M4y,, 8/4 F_4/o(7,)|-¥2 c.g.8. units, . (8) 


where «,,=m,,/mMo, Mo being the actual electronic mass. 

It suffices to insert for 7,, in (8) the ordinary average value ignoring 
the space-dependence of the electron density. This is determined by 
the relation 

n= 2h*(2rm kT) PF 1)9(n); 
Ke) F 1j(,)=2'070 x 10-7 8/20, 8/2 ¢.9.8. units. en O) 

Extensive tables of the Fermi—Dirac function ¥ appear in Dingle 
and Roy (1955). In connection with the present paper we reproduce 
here a short table of F4/. and F_4)5. 


4) F 120) F -1/2(n) 7 F 127) F ~4)2(n) 
—4:-0 0-0182 0:0181 3°5 5-458 2-026 
—3°5 0:0299 0:0296 4:0 6-512 2-186 
—3:-0 0:0489 0:0481 4-5 7:642 2-334. 
—2°5 0:0798 0:0776 5:0 8-844 2-473 
—2:0 0:1293 0-1237 5:5 10-11 2-604 
—1-5 0-2074 0-1933 6:0 11-45 Dien 
—]-0 0:3278 0:2940 6:5 12-84 2-845 
—0:5 0:5075 0:4312 7:0 14:29 2:957 

0:0 0:7651 0:6049 1:5 15:80 3:065 

0:5 ifoilalyl 0:8077 8-0 17:36 3°169 

1-0 1-576 1-027 8:5 18-96 3:270 

1-5 2:145 1-249 9-0 20-62 3:°367 

2-0 2-824 1-464 9-5 22e30 3:461 

2°5 3°607 1-666 10-0 24-08 3°553 

3:0 4-488 1-853 
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For highly degenerate electron systems, 7,>0, application of the 


asymptotic relation 
F io) SE) 


leads to a value for the screening distance 


Feaegenerate= (77/3)¥/8h y/x/4rrev/m,, Shoe 
= 3-665 X 10-5«1/2n-M6y 1/2 o.g.s. units, . . (11) 


in agreement with Mott (1936, see also Mott and Jones, p. 87, where 
their q~! equals our &). 
For non-degenerate electron systems, 7,, <0, the approximation 


FA ae. id tee ee ee 
leads to the expression 


R (xkT'/40n)/?2/e=6-902K1/27'1/2n-1/2 6.9.8. units. . . (13) 


We now consider two typical examples. Suppose a specimen of 
silicon contains some phosphorus atoms as impurities. Phosphorus 
has a valency of five and silicon of only four; in order to permit the 
phosphorus bonds to fit in with the fundamental structure of the silicon. 
based on four covalent bonds, one bond of each impurity phosphorus 
atom is destroyed. The result is that each phosphorus atom acts as 
a single positive charge and—at least at reasonably high temperatures— 
provides one electron to wander through the crystal. The degeneracy 
temperature of the resultant electron system is given by the customary 
formula (e.g. Mott and Jones 1936, pp. 54, 180) 


non-degenerate 


T —s he SAG 2/3 — 4.93] x 10-112/8y-1 i 
degen. = Bf (<) nei = 4 231% n**a-* ¢.g.8. units, . (14) 
where » is the number of free electrons per unit volume and «=m/imp, 
as before. The number of silicon atoms per cm? is 5-00 x 1022. 

For our first example, we suppose that one impurity P atom is 
introduced for every 100 Si atoms, so that n=5-00x 102°em-3. If we 
assume that the electrons behave as though perfectly free, so that «=1, 
then 7'jogen—= 2665°K, showing that with this large impurity content the 
electron system is totally degenerate at room temperature. Since for Si 
k=11-9 (Briggs 1950), we deduce from eqn. (11) that the screening 
distance R is 4:488A. For comparison, the cut-off distance implicitly 
assumed in the theory of Conwell and Weisskopf (1946, 1950) would be 
3n-1/3—6-300A. We should, however, remark at this point that the 
validity of the theory is somewhat dubious for such high impurity 
concentrations as that considered in this example. For instance, the 
screening distance R=4-488 A turns out to be not quite twice the average 
spacing of the Si atoms, N,g~/%=2-714A, so that the insertion of the 
bulk value of the dielectric constant « of Si is questionable. It is also 
possible that quantum exchange and correlation effects, neglected in 
the present paper, may influence the results appreciably. 
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For our second example, we suppose that only one impurity P atom 
is introduced for every 10° Si atoms, so that n=5-00x 1017 em-3, Then 
T degen. = 26-65°K, showing that the electron system is now quite non- 
degenerate at room temperature. Equation (13) then gives, at a 
temperature of 17°c, R=57-344. For comparison, the Conwell—Weisskopf 
cut-off distance would be 63-00 A. 

According to our formulae the screening distance is not simply pro- 
portional to the spacing of the impurity ions N,-1/ even in the 
limiting cases of strong and weak degeneracy, so that customary cut-off 
procedures which assume this are strictly speaking unreliable. Their 
apparent success is clearly explained by the largely fortuitous rough 
numerical agreement under normal experimental conditions between the 
values of the theoretical screening distance R and the assumed cut-off 
distance 4N,-1/. 

Screening of Acceptors by Holes 

Writing the statistical distribution function for electrons of energy ZL, 

measured upwards from the bottom of the appropriate band, as 


: (15) 


fe exp (BE PTF 
that for holes (=missing electrons) must be 
1 
fee expC, Eu EET eae 


where #,, is the energy the missing electron would have had. Since elect- 
rons will normally be missing only from the top of a band, it is convenient 
to measure the energy LH, of holes downwards from the top of the appropri- 
ate band by writing ree sre ee ee 


where A is the energy gap between the top of the band supplying the 
holes and the bottom of the band supplying the electrons. In other 
words, —A is the energy at the top of the hole band as measured on 
the energy scale used for the electrons. We may then write, in complete 


analogy with (15), 1 


Jo= exp (Hy — bp) RT TF 1’ oe 


where €,——A—{,,. 

‘Adopting this notation, itis clear that the formulae giving the screening 
distance when acceptors are surrounded by holes are of just the same 
form as those for the screening distance when donors are surrounded by 


electrons. Thus by analogy with (7), (8) and (9), 

b(r)=— —exp (Ear Rye. 3 wanerea oye (19) 

1677?e?m,,3/2(QrkT)\/? 
Kh3 


and 7, is determined by the relation 
p=2h-3(20m,kT PF yjolqy), -~- + - +» (21) 


pe R= Ba ee Fe 2)"(20) 
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» being the effective hole mass, p the average number of holes per 
Ant volume, and ,={,/kT7=—A/kT—n,. Limiting formulae for 
complete degeneracy and non-degeneracy may similarly be written 
down at once by analogy with (11) and (13). 


Screening of Donors (and Acceptors) by both Electrons and Holes 


It will be seen from (5) that the space-dependence p(r) of the screening 
charge round a positively charged donor ion is proportional to e?, and 
to this extent is therefore independent of whether these charges are 
electrons or holes. The physical explanation of the screening of donor 
ions by holes of like charge is that since these holes are repelled from the 
donor, any incident particle must approach from a region which is already 
filled with positive charge, and is therefore less affected by the fixed 
positive charge of the donor, only changes in potential energy being ~ 
effective in scattering processes. 

For the screening of donors and acceptors by electrons and holes, we 
have 


é e 
Paonor= pe peedy (—r/R) ; Pacceptor= — eae (—1/f), . (22) 
where for either donors or acceptors 
16772e2(2rkT')1/2 geek 
Be [22 Fale tmp? Friel) (23) 


§3. EVALUATION OF THE FERMI-LEVEL 


The fundamental quantity »,—€,/k7' (the subsidiary quantity », 
being just given in terms of this by the relation »,——A/kT'—n,,) may 
be calculated from the condition of overall electrical neutrality 
a—Vap ee ee ELS 
in which » and p are given in terms of 7», (or »,) by (9) and (21) 
respectively, and 

Na a 

ls 2 exp [n, +H /kT]+1’ Nate, 2 exp [y,+H,/kT]+1° * ) 
Here Ng is the total number of donor atoms introduced, comprising N 4; 
positively charged and Na—N4, neutralized by bound electrons ; JN, is 
the total number of acceptor atoms introduced, comprising V,; negatively 
charged and N,—WN,; neutralized by bound holes ; Hy is the energy gap 
between the bottom of the electron band and the donor levels, namely 
minus the energy of the donor levels as measured on the electron energy 
scale, i.e. counting upwards from the bottom of the electron band reckoned 
as of zero energy; and —#, is the energy of the acceptor levels as measured 
on the hole energy scale, i.e. counting downwards from the top of the hole 
band reckoned as of zero snares Note that #, (hole scale)=A—E, 
(electron scale). The factor 2 multiplying the exponentials i in (25) arises 
from the fact that bound electrons and holes may be held in either of 


n—p=N 
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two spin orientations—see, for e.g. Shockley (1950, pp. 248, 475), Lands- 
berg (1952, 1953), Guggenheim (1953). 


$4. Tor Time or RELAXATION 


In all cases considered in the present paper, the magnitude of the 
potential energy between an incident electron or hole and a charged 
ionized donor or acceptor may with sufficient accuracy be written in the 
form 


2 
| Vir) | = |ed(r) |=<exp(—7/R). .. (26) 


The resultant differential elastic cross section J(#), defined as the number 
of particles per unit time scattered through an angle @ into unit solid 
angle when on the average one particle is incident with velocity v on 
unit area per unit time, is given by Born’s approximation (e.g. Mott and 
Massey 1933, p. 119): 


82m (@ 
V1(0)= Fa | 


where K=4amv (sin $0)/h. 
An elementary integration then yields 


sin Kr | V(r) | r dr ies ha eR 
0 


e2 

Ws cmv" {2 sin? 40+ (h2/87?m?v)/R?} 

As is well-known (e.g. Mott and Jones 1936, p. 262), the reciprocal 7~1 

of the time of relaxation is equal to the total probability per unit time 

that a particle is scattered, except that a factor (1—cos@) must be 

introduced (as in the theory of gaseous diffusion) to give extra weight 

to large-angle collisions. Thus if there are N, charged ions per unit 
volume, en 

71 Nw 20 | (1—cos 0)1(6) sin 0 dd. Pea (29) 


0 


(28) 


Taking (1—cos 0)=z, this reduces to 


hg a a Oh 
een 0 | o [e-+ (h?/87r2mm2v2)/ R?]? 
27N,e* ie 
=e [mate]. Seren) 


where =1672mv?h2/h2, in agreement with Mott (1936) provided we 
replace his Z? by x-*. Written in terms of the energy H=}mv* of the 


incident particle, ean od . 
rt inane), eee 2" (32) 

where g(€)=In (1+ €)—€/(1+€) ae) (33) 

with E=32r?m R7E/h?. Meet set (ae) 


The form of these results, and of R, is independent of whether the 
seattering centre is a donor or acceptor ion. For incident electrons 
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m must be understood to refer to m,, while for incident holes m=m,. 
The abbreviation In refers to the natural logarithm. 

It is convenient to define the ‘ temperature’ of an incident particle 
by the relation H=kT partic: noting that 1ev=11 605°%. Then in 
¢.g.s. units, 

E=9-049 x 10140? 7 particle 5%. Cee WO) 


where «=m/m,. For incident particles of very low energy, €<l, 
g(é)=4€2, and hence 7-1 cc \/E# as for scattering by thermal agitation. 


§5. APPROXIMATION TO THE TIME OF RELAXATION 


Owing to the comparative complexity of the function g(£) defined 
by (33), it is necessary to replace it by an approximation before 
performing the required integration over energies, except of course 
when the electron system is sufficiently degenerate for the Sommerfeld 
(1928) asymptotic series to be applicable. Let us first consider the 
likely magnitude of € in practical cases. 

In the first example considered in §3 (one P impurity atom for 
every 100 Si atoms), the electron system was degenerate, so that 
T particle /'aegen. = 2665°K. Inserting the value of the screening distance 
R=4-488 x 10-8 cm previously determined in §3, it is found that the 
only important values of € lie near €)=4-857. 

In the second example (one P for every 10° Si), the electron system 
was non-degenerate. The incident particles-must then be travelling with 
thermal velocities, so that Tyartice~Z'. For reasons which will become 
clear later in this section, we shall take 7',,tice=3°6727' when con- 
sidering the most important problem, that of the electrical conductivity. 
Putting 7=17°c=290°K and R=57-:34x10-8§cm as deduced in §3, 
the important values of € lie around €,=316:8. 

As already remarked, for degenerate samples the Sommerfeld (1928) 
asymptotic series is applicable, and the complete function g(é) may be 
used without undue complication. Hence the conclusion to be drawn 
from these figures is that in all normal circumstances in which an 
approximation to g(€) is really necessary, € is considerably greater than 
unity for incident particles with average energy. Under such conditions, 
g(€) is such a slowly varying function of €, and therefore of HZ, that it 
suffices to replace it by a constant of magnitude g(&). 

In order to assess the error introduced by such an approximation, 
we consider the formula for the electrical conductivity of a non-degenerate 
electron gas. This involves the integral 

+00 2) 
B87 exp(—E/kT) dE x | a exp (—x) g-(Eou/3-672) dx, . (35) 
J0 “0 
where #=T7'yarticte/Z'- The problem is then to estimate the error caused 
by taking the slowly varying quantity g~1(&)2/3-672) out from under 
the integral sign and inserting some suitable value of 2 in its argument. 
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Now it is easily found that a quarter of the total value of 


| x? exp (—x) dx 
0 
is contributed by values of «<2x,)4=2-537, a half by C0, 3 OTe 
(and hence the introduction of this figure—the median—earlier in the 
section), and three-quarters by x<a3,—5:101; for these particular 
x-coordinates the factor g~1(£,x/3-672), with €,—316-8 as in our second 
example, takes the respective values 0-2274, g~1(&))=0-2099, and 0-1964. 
The error committed by assuming that the factor g-! appearing in the 
integral may be replaced by g~1(&,), its value at the median (the bisector 
of the area under the curve of the integrand) will be of the same order 
of magnitude as the difference between the results of such an assumption, 
and of those obtained by replacing g-! by g~'(£o4/4/%1/2) over one-half 
of the range and by g~1(£9%3,4/%1/2) over the other. Thus the error will 
be of the order of 


0-2099— 4(0-2274-+ 0-1964) 


Oe eo 
0-2099 > 100%, By 


We believe that the insertion of the particular value x=2 qian IN gut: 
where 2.aian=1/2 represents the median of the integral remaining 
after g~1 1s taken outside the integral sign, will under most circumstances 
give very much more reliable results than the customary insertion of 
Lmax, the coordinate for which the remaining part of the integrand 
possesses a maximum (e.g. Conwell and Weisskopf 1946, 1950, Shockley 
1950, p. 279), since our procedure automatically compensates for the 
asymmetry of the remaining integrand. Indeed, our method can still 
be used even if the remaining integrand possesses no maximum, as 
sometimes happens (Dingle and Roy 1955). 


§6. COMPARISON WITH FORMULAE OBTAINED BY CuT-OFF PROCEDURES 


In view of the frequent use of cut-off approximations for interpreting 
data in the literature, it is particularly interesting to compare our results 
with those obtained by such means, taking as a basis for comparison 
the formulae quoted by Shockley (1950). The quantity 7N ,e4H~*/?/(2m)'/?x? 
appears as a multiplying factor in both our expression and that given by 
Shockley* (1950, p. 264). For a semiconductor containing charged donor 
ions and electrons, we have therefore to compare our equation 


g(€)=In (1+ €)—E/(1+€) 
where € takes the limiting forms 


gr 1/8 SrKmkT 
ens 


er) = ee 
Casrencrate =2 ( 3 E and ron degenerate e2 h2n 


Ce ta a 
* Wilson (1953, p. 266) has an additional factor $ in 77?. 
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with Shockley’s equation 4(&)=In (1+ &) 


where* é K E. 


P 


: e2n 1/3 
It is at once seen that in the limit of complete degeneracy 
&/E=}(3/7)3=0-4924, ) 
a pure number. Our formulae then agree with Shockley’s at least as far 
as to which variables are involved. For our first example, with one P atom 
for every 100 Si atoms, we found in § 5 that €)=4-857, giving @)=2-392, 
g( €o)=0-9384 and 4(&))=1-905. We conclude, therefore, that for typical 
degenerate electron systems Shockley’s formula gives values of 7-1, and 
therefore of the resistivity, too high by a factor of about 2. We may 
also note at this point that owing to differences in their cut-off approxi- 
mation for the degenerate case, Johnson and Lark-Horovitz (1947) 
obtain a value of 7! smaller than Shockley’s by a factor 7: their 
resistivities would therefore be too low by a factor of about 3/2. 
In the non-degenerate limit, 
6 /E=h?n?!3/8amkT =1-389 x 10—-140—-1n?47—! in c.g.s. units, 

and any agreement between the results of our theory and those of 
arguments based on. the customary cut-off approximations must be 
largely fortuitous. For our second example of one P atom for every 
10° Si atoms, &/€=0-03018 at room temperature. Thus while €)=316-8 
as found in §5, &)=9-560. These correspond to g(&))=4:765 and 
4( 6 )=4:526, so that although cut-off procedures for non-degenerate 
electron systems appear to lead to formulae not even involving the 
correct powers of the variables, for the typical experimental conditions 
discussed the numerical results obtained by their use fortunately turn 
out to be reasonably accurate. 
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ABSTRACT 


By combining results obtained from a methane counter with those 
from nuclear emulsions the cross-sections for the !2C(y3«) reaction at 
17-6, 14-8 and 12-3 Mev are found to be (1-70-40-24), (0-33-+.0-07) and 
(1:15+-0-6) x 10-78 cm? respectively. Analysis of the data suggests 
a strong resonance in the excitation function near 12-3 Mev and possible 
resonances near 15 and 16 Mev. 


$1. InTRODUCTION 


CONSIDERABLE attention has been given to the reaction !C(y3«) and 
the position has been reviewed recently by Titterton (1955). Many 
observers have determined the cross section for the “Li(py) resonance 
radiation by the photographic plate method, a knowledge of this being 
important since it is often convenient to determine cross sections for 
other, rarer, reactions occurring in the emulsion in terms of it. The wide 
spread of values which have been reported (table 1 and fig. 3) is un- 
satisfactory and makes a redetermination of the cross section at 17-6 Mev 
desirable ; at the same time it was hoped that by using other lines in the 
7Li(py) spectrum, particularly the broad one at 14-8 Mev, information 
could be obtained on the resonance absorption of radiation by states 
of 12C, evidence for which had been presented by Goward and Wilkins 
(1953). 

To determine the cross section at 17-6 Mev a methane filled proportional 
counter was used since this enables data to be collected more quickly 
and with better energy resolution than is possible by the emulsion method. 
However, the information obtained is less definite; only the total 
disintegration energy is measured and it is impossible to distinguish 
between events where one particle strikes the counter walls and events 
due to y-rays of energy lower than the main line, although appropriate 
corrections can be applied. For this reason the proportional counter 
was used to obtain the absolute cross section at 17-6 Mev while the 
photographic plate method was employed to investigate the reaction 
at lower energies and to determine cross section values at 14-8 and 
12-3 Mev relative to the 17-6 Mev measurement. 


* Communicated by the Authors. 
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§ 2. DescRIPTION AND PERFORMANCE OF COUNTER 

The counter is cylindrical, of length 30cm and internal diameter 
9-1em. A tungsten wire of diameter 40 is taken through a Kovar-glass 
seal inside a 1-5 mm diameter copper capillary tube, hangs along the axis 
of the counter and is kept taut by a nickel weight ; 18 cm of wire is 
exposed within the counter. The outer pressure chamber is made from 
stainless steel 1/16 in. thick and lined with a graphite cylinder and end- 
plates, also of 1/16 in. thickness, which form the cathode surface of the 
counter. This lining prevents charged particles, produced in the steel 
walls, from entering the sensitive region. The only important reaction 
occurring in the graphite at these energies is the reaction under study and 
consideration of the maximum depth at which a photodisintegration 
event could occur and still lead to the dissipation of more than 5 Mev 
in the active volume shows that <0-1°%, of observed pulses are due to 
reactions occurring in the graphite walls. Attached to the main chamber 
is a gas purifier of the circulation type containing evaporated sodium 
metal. This reacts strongly with water vapour and oxygen, the two 
most objectionable of the common electronegative counter gas con- 
taminants, but does not react with methane. Gas is pumped in and out 
of the counter through two needle valves, the moving parts of which are 
enclosed in bronze sylphon bellows to prevent grease from coming into 
contact with the gas. 

The fillimg procedure was as follows: the air in the counter was 
replaced by methane and a pellet of sodium dropped into the purifier. 
The methane was then pumped out until a pressure of ~10-4mm Hg 
was obtained when the sodium was evaporated on to the walls of the 
purifier tube. Pumping was then continued for half a day, the chamber 
being outgassed periodically by heating with a flame. Methane was then 
allowed to pass slowly through a copper coil immersed in solid carbon 
dioxide and acetone to freeze out water vapour and was admitted to the 
counter until the pressure reached 2:38 atmospheres. 

At this pressure and at the working voltage of 2-4 kv the gas multiplica- 
tion of the counter was found to be 3 times. The counter pulses observed 
on an oscilloscope had rise times between 10 and 100 pusec—a rather 
greater spread than expected. According to measurements by English 
and Hanna (1953) of electron mobilities in methane, the spread due 
to electron straggling would have been expected to be ~3 psec; this 
difference was not investigated and is not understood. To handle the 
long pulses satisfactorily rise and clipping time constants of 80 sec were 
used in the amplifier, the output of which was fed to a 120 channel pulse 
amplitude analyser. 

An energy calibration for the counter was obtained using the spectrum 
of «-particles from Th-X and its daughter products decaying with a 
half-life of 3-65 days. A small amount of polonium was included in the 
counter to give a continuous check on the calibration. The energy scale 
was found to be linear to within 1%. 
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§ 3. THe CouNTER EXPERIMENT 


Resonance radiation from the 7Li(py) reaction was excited by 
bombarding a thick lithium metal target with 500 kv protons, the counter 
being situated with its axis perpendicular to the proton beam, and 
centre 10-8 cm from the target at 0°. Irradiation was carried out in 
8 hour runs, alternating with background runs of similar duration ; 
the C(y3x) counting rate was small (~140 per hour) necessitating a 
careful measurement of the counter background. Actually the C(y3«) 
peak produced by the 17-6 Mev y-rays yields «-particles with a total 
energy of 10-3 Mev, well above the energy of the background due to 
radioactive contamination, and in 44 hours no background pulses over 
9:2 Mev were observed. However, a knowledge of the background, 
mainly from polonium (100 counts/hour) and radium (20 counts/hour), 
was important in an analysis of the tail of the 17-6 Mev peak and of peaks 
due to lower energy y-rays. 

The radiation was monitored using a thick-walled brass Geiger-counter 
identical in construction with that described and calibrated by Barnes 
etal. (1952). Readings of pulse counts were recorded at one hour intervals. 
The relative photodisintegration yield per quantum in the four runs was 
POOR I-00) 1206750206. 


§ 4. THE Cross SECTION AT 17:6 Mev 


The pulse spectrum resulting from 31 hours irradiation is shown in 
fig. 1 (a) which, after subtraction of the background given in | (0), leads 
to the curve of 1(c). The high energy peak contains approximately 
2900 counts, has a half-width of 0-8 Mev, and centres at 10-30--0-06 Mev 
as expected from the mass data which gives the reaction threshold as 
7-28 mev. Contributions to the spectrum below the main peak come 
from y-ray pile-up pulses, «-particles from 17-6 Mev events which are not 
contained entirely within the counting region, and '7C(y3x) events induced 
by y-rays of energies less than 17-6 Mev. The contribution from y-ray 
pile-up pulses was estimated by extrapolating the low energy tail measured 
with the kicksorter bias at a reduced value. The ‘ pile-up’ pulse contri- 
bution, although small was significant, amounting to 130 counts per 
channel at 5 Mev, decreasing to 10 counts per channel at 8 Mev for a 
31 hour run. 

To allow for losses to the walls a knowledge of the ranges and directions 
of the «-particles is necessary. The correction applied assumes the 
disintegrations to be two-stage : 


12C_>%+ 8Be ; 8Be+>2a 


and to be isotropic, with 93% of the transitions through the 2-9 Mev 
level of 8Be, and the remaining 7° through the ground state. A range— 
energy relation for «-particles in methane was deduced from the curves 
of Livingstone and Bethe (1937) and correction factors were derived from 
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those for single «-particles to give the fraction of events in which one or 
two «-particles of a star intercept the boundaries of the counter. The 
energy distribution of the resulting degraded pulses was calculated and 
amounts to approximately 42 counts per channel between 6 and 9-5 Mey. 
The shape of the distribution has little effect on the magnitude of the peak 
at 17-6 Mev, but makes the yield at lower energies uncertain and it is for 


Fig. 1 


c 
eee Renne! 


Counter experiment. 
(2) Pulse spectrum obtained in 31 hours irradiation of CH, counter. (b) Counter 
background spectrum normalized to 31 hours running time. (c) Pulse 


spectrum resulting from disintegration of 12C (difference of (a) and (b)). 
(d) Corrected distribution. 
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this reason that the cross section at 17-6 Mev only has been calculated from 
the counter data. 

The final distribution fig. 1 (d) was derived from fig. 1 (c) by subtracting 
the * pile-up’ pulse distribution and the pulse distribution arising from: 
wall losses and then adding the loss corrections at the appropriate energy 
to the 17-6 Mev peak. The peak is asymmetrical probably because the 
amplifier time constants were shorter than the rise-times of some of the 
pulses ; the value 80 sec was chosen for the time constants to give the 
best compromise between lack of uniformity in pulse height and difficulties 
due to the y-ray ‘ pile-up ’ pulses. 

A small correction for .absorption of radiation in the counter walls. 
was calculated from data given by Heitler (1949). The intensities of the 
y-ray lines were assumed to be in the ratio 17-6: 14:8: 12:3:: 0-62+5% : 
0:36+12%: 0-02+50%, the first two being based on data of Stearns. 
and McDaniel (1951) the last figure being based on Inall (1954) and 
Titterton (1953). 

The final value for the cross section at 17-6 Mev is 


F17.6(127Cy3a)=(1-70-L 0-24) x< 10-28 cm2. 


Table 1. Absolute measurements of the cross section for ?C(y3«) 


Author Gmean Cine 
Waffler and Younis (1949) 0-8 +0-3 x 10-28 em? 
Glattli, Seippel and Stoll (1952) | 1-75+0-25 2-4 10-78 cm? 
Goward and Wilkins (1953)* 1-4 +0-26 1-70 40°34 
Greenberg, Taylor and Haslam 
(1954)+ diied 
Present work 1-19 +0-20 1:70 40-24 


* Combination of lithium y-rays and bremsstrahlung measurements. 

+ Bremsstrahlung measurement. yaa 

The value for o,7., from Li y-rays depends on the relative intensity assumed. 
for the 17-6 Mev line which is not the same in the experiments listed. 


§ 5. PHoTOGRAPHIC PLATE MEASUREMENTS 


Incidental to another investigation some 1000 C(y3x) ‘stars’ were 
measured in Ilford 200 » type E, boron-loaded emulsions which had been 
exposed to 7Li(py) resonance radiation. The stars were identified by the 
usual method of momentum balance and the number-energy histogram 
is given in fig. 2 (a). Three peaks are resolved, the main one at 17-6 Mev 
contains 742 events, the broad 14-8 Mev group has 86 events while the 
group centring on 12-5 Mev has 17 events. These numbers may be 

compared with 2160, 300, 35 observed by Nabholz et al. (1952) and 390, 
50, 12 observed by Goward and Wilkins (1953). 
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Loss corrections were made by the geometric method of Goward and 
Wilkins (1953) and the smooth distribution (full line of fig. 2 (a)) obtained. 
The final result is compared with the methane chamber data in fig. 2 (0), 
the curves being normalized by the area of the 17-6 Mev peaks. 


Fig. 2 


N/ 25 MeV 


Units 


Arbitrary 


Photographic plate. 

(a) Histogram of energy released in 845 #C(y3q) ‘ stars’ with corrected curve 
superimposed. (6) Comparison of results from methane counter and 
photographic plate normalized by areas of the high energy peak. 
* counter, o plate. 
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§ 6. Cross SECTIONS aT 14:8 AND 12:3 Mev 


Taking the cross section at 17-6 Mev to be 1-70 10-28 cm? and the 
relative intensities of the 17-6 and 14-8 Mev lines to be as given in § 4 the 
photographic plate data gives : 

6 14.8(12Cy3a)=(0-33-40-07) x 10-28 em2. 

The group centring on 12-5 Mev in fig. 2 (a) could be due to a reaction 
of the type #C(yy’)3«, if a level of 1C is available at the appropriate 
energy. However there is strong evidence from experiments in this 
laboratory (Titterton 1953, Inall and Boyle 1953, Inall 1954) that there 
is a weak line at 12-3 Mev corresponding to a transition to a state of 
8Be at 5:3 Mev. We therefore discard the possibility that the group may 
be due to the inelastic scattering of y-rays and calculate the cross section 
for the C(y3«) reaction at 12-3 Mev taking the intensity of these y-rays 
to be (2+1)% of the total radiation. The figure obtained is : 

O753(2 -Gysx)—(1-'15420-6) x 10- 2° em?. 
Using the same set of conditions and normalizing all cross sections to our 
17-6 Mev figure, cross sections at 14-8 and 12-3 Mev have been calculated 
from the histograms published by Nabholz e¢ al. (1952) and Goward and 
Wilkins (1953) and are presented in table 2. It will be seen that there is 
good agreement among the three experiments for the cross section at 
14-8 mev but considerable differences at 12-3 Mev. 


Table 2. Cross section for #2C(y3x) at 14:8 and 12-3 Mev, normalized to 
O17.¢= 1:70 10-28. The lower two sets of results were calculated 
from the published histograms. The relative y-ray intensities 
used are : 

1(17-6 Mev) =0-6245% (14-8 Mev) =0-36-+12% (12:3 Mev) =0-02-50% 

nnn ee SS 
| 


K etl iv : . ‘ 
= Cross section x relative Cross section in 10~?8 em? M. 
Source y-ray intensity ean 


17-6 Mey | 14:8 Mev 12-3 Mev 17-6 Mev 14:8 Mev | 12:3 Mev 


Present work, 


CH, counter | 1-05 1:70-0:24 
Present work 1:05+2-3]0-118+7 |0-023+16 0-33-40-07 | 1:15-L0-6 | 1-19-+.0-20 
Nabholz et al. 1-:05-+1-4 | 0:144-+3-7| 0-015+11-5 0:40-+0:8 | 0-75 1-21+0-20 


Goward and : 
Wilkins 1:05+3-2]0-130+9 |0-032+20 0:36-L0-07 | 1-6 1-21+0:-20 


_ eee eee ee ne 
The high value of o 5.3 implies the presence of a !2C level in this region 
and the resonance absorption of radiation of this energy. Such an 


interpretation is supported by some of the bremsstrahlung data and by 
results obtained in a study of the reaction #C(nn’)3« by Livesey and 
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Smith (1953). Excitation functions over the region 12-18 Mev have 
been measured with poor agreement by five groups and are presented 
in fig. 3; the present three points have been plotted for comparison. 

Inspection of both the methane chamber and the photoplate data 
over the width of the 14-8 mev line, fig. 2 (b), shows elevation of the points 
near 15 and 16 mev and a trough at 15-5 Mev as indicated by the arrows. 
Taken separately little weight could be placed on either result, but the 
two sets of data support each other. Moreover, examination of the 
photoplate results of Goward and Wilkins (1953) and Nabholz et al. (1952) 
show the same features to be present in each of these experiments. 
Taken together these results lend support to the presence of peaks in the 
cross section at 15 Mev and 16 Mev as observed by Goward and Wilkins 
(fig. 3) but suggest that the resonances are more marked than their 
excitation function indicates. 


Fig. 3 


Goward € Wilkins 
(1953) 


Dowson ¢ Bigham 
(1953) 
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Millor ¢ Cameron 
(1953) 


Ed Telegdi 
mh hear 


Cross- Sectio 
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2 14 16 18 EY.Mev 


Cross sections for the reaction 12C(y3«) over the region 12-18 Mev. #% present 
results plotted. 


The present experiment therefore gives strong evidence for the resonance 
absorption of y-radiation by the !2C nucleus near to 12-3 Mev and suggests 
the presence of other resonances near 15 and 16 Mey. 
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SUMMARY 


Thirring’s calculation of the gravitational field near the centre of a 
rotating spherical shell is corrected and elaborated. First, an additional 
term, representing elastic stresses, is introduced into the energy-momentum 
tensor; without this term the conservation laws cannot be fulfilled and 
the shell could not keep together during rotation. Secondly, the 
dependence of the gravitational effects on the mass-distribution over the 
shell is investigated. The quasi-centrifugal forces are to be attributed 
entirely to the mass-distribution on the shell arising from the mass- 
increase with velocity. 


§ 1. INTRODUCTION 


In a classical paper, H. Thirring (1918, with minor corrections 1921) 
calculated the gravitational field near the centre of a rotating thin 
spherical shell, using Einstein’s field equations in the linear approximation 
(Einstein 1916). Thirring found that the rotating shell produces near 
its centre forces analogous to the Coriolis and centrifugal forces which 
arise in rotating reference frames in classical mechanics—-his aim was 
to confirm by an example that Einstein’s theory embodies Mach’s 
principle of the relativity of inertia. 

However, as Professor C. Lanczos has kindly pointed out to us, 
there is an inconsistency between Thirring’s assumptions. Thirring 
assumes that the material of the shell is incoherent, so that its energy- 
momentum tensor has the formt 

Tix pulyt as (1) 
where p is the density, v‘=da"/ds the 4-velocity of the matter composing 
the shell. Thirring proceeds to calculate approximately the gravitational 
field near the centre of the shell, working to the second order in the 
angular velocity w of the shell, and neglecting the self-gravitation of 
the shell itself. But to this approximation, the conservation law$ 


TOY ose ¢ 
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* Communicated by the Authors. 

+ Private communication. 

{ Greek suffixes range and sum over 1, 2, 3, 4. The units are grams, centi- 
metres and light-seconds, so that c=1. The coordinate x4=it is purely 
imaginary, the other coordinates real, and the metric has signature +4. 

§ A semi-colon (;) denotes covariant differentiation. 
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is not satisfied ; if (2) is to be satisfied by (1), the incoherent particles 
cannot move in circular paths, but must follow geodesics, which is to 
say that a shell of incoherent matter, if set in rotation, will fly apart 
along the tangents to these paths. To say that the shell should be 
regarded as composed of distant stars or galaxies is not to resolve the 
difficulty, because if there is no interaction between them then such 
stars or galaxies cannot follow the required paths. 

This paper is concerned not with the cosmological significance of the 
problem, but only with the resolution of the inconsistency. This will be 
_ effected by supplementing the energy-momentum tensor (1) of incoherent 
matter by a term representing the ‘elastic’ interaction between particles 
of the shell. 

In general, such a substitution would require a generally relativistic 
elasticity theory, but in the present case the symmetry and simplicity 
of the problem allow us, under natural physical assumptions, to 
determine uniquely the contribution of such an elastic interaction to 
the energy-momentum tensor, without developing such a theory. 


§ 2. DEFINITION OF THE ENERGY-MOMENTUM TENSOR 


We thus adopt for the total energy-momentum tensor of the material 
of the shell ) 

T= pv'y’ + BY, SAL eeye release {3) 
where the added term EH’ arises from elastic interactions. 

It is clear from the discussion in § 1 that what is needed to hold the 
shell together is just a tangential stress along the parallels of latitude. 
One could simply assume this from the outset, and determine the 
magnitude of this stress from the conservation laws, thus arriving 
directly at eqns. (14) and (20) below. Alternatively, one may proceed 
from general assumptions about the properties of H"” (which would 
apply also in other circumstances), and show that under these assumptions 
only the tangential stress term affects the final results, in the approximation 
considered. 

This alternative procedure will now be carried out, but the essential 
results (beginning with the line below eqn. (19)) can be understood 
without going through the details of the calculation. 

The elastic term E“’ is assumed to have the following properties, 
which are natural generalizations from classical elasticity theory : 

(i) It is symmetric : 

BY — Bh, fe Se ay (4) 


This is perhaps scarcely necessary to mention. 
(ii) It has no component normal to the shell : 
Ben ,=0, eee = ar (5) 


where at any point of the shell , is a 4-normal to the world tube (3-surface) 
of the shell. This generalization of the corresponding condition in the 
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ordinary theory may be regarded as an application to a thin shell of 
the continuity conditions of Synge and O’Brien (1952). Its physical 
meaning is that there is no transfer of elastic energy-momentum across 
the surface of the shell. Since v lies in the world tube 3-surface, 
v'n,—=0, so that 

Tens, 
as one would expect. 


(iii) It has no component in the direction of motion of the shell : 
EYv,=0. othe ea tgs te ee 


There will in fact be a contribution to the energy-momentum tensor 
arising from the elastic energy. By considering an element of the shell 
in a local reference frame in which it is instantaneously at rest, one can 
see that this contribution is of the form ev“v’, where ¢ is the elastic energy 
density. This contribution may be incorporated into the kinetic 
energy-momentum term pv“v’, and the remainder will naturally satisfy 
eqn. (6). In general circumstances the term ev“v” could not be neglected 
—the material density would have to be increased to take account of it— 
but in the considerations of § 3 it will be neglected, because it is pro- 
portional to the square of the stress and therefore of order w*, where 
w is the angular velocity of the shell, while the calculations are to be 
carried out only to order w?. 


(iv) Its physical components, in a local reference frame in which the 

relevant point of the shell is instantaneously at rest, are bounded : 

Eh, Oh, O <A, ate okey See tae 
(A a constant for the whole shell), where 4,‘ denote v, and any three 
unit vectors orthogonal to v, and to each other. 

Under these assumptions, #” can be determined uniquely from the 
symmetry and steady state of the problem. This is a little surprising, 
since nothing has been assumed about the elastic properties of the shell. 

Following Thirring, we consider a thin spherical shell of radius a, 
rotating with constant angular velocity w in flat space-time (deviations 
from flatness due to the presence of the shell are neglected in the 


neighbourhood of the shell itself). The metric in spherical polar 
coordinates* is 


ds*= (da*)?- (dr)?-+-712(d0)?-+-r? sin? 0(dd)?. PE gL Ai a 


In this coordinate system, the velocity 4-vector of any point of the shell, 
for a uniform rotation about the axis 0=0, is given by 


oy sora ve eos e: _ dg w 
ot = dat 9 vt dat” ip ndacteen ths sie! 
so that 
Y=—iwv', (vt)?—wa? sin? 6(v3)2= — 1, Se LO) 


* Labelled 1, 2, 3, 4 in the order 1, 6, 4, x*. 
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and ‘so 
7p), v=wy(9), vt=iy(6) | 
; ek en 
x(0)= (1 —w?a? sin? 9)-1/2, i 
In this coordinate system, the unit normal n, at any point of the shell 
thas the covariant components (1, 0, 0,0). Unit vectors A» Hy, in the 
4 and ¢ directions have respectively the covariant components (0, 7~1, 0, 0), 
(0, 0, 7-1 cosec @,.0). Since any point of the shell can be brought 
instantaneously to rest by Lorentz transformation to a frame with 
relative velocity wa sin 6, it follows from (7) that the physical components 
of H’” must be bounded also in the reference frame defined by (8) (that 
is, in the Minkowskian frame in which the shell rotates about a fixed 
axis). 

We now write down the conservation eqns. (2), simplified with the 
aid of conditions (i)—(iv). In these calculations, the density p (which 
is constant in Thirring’s work) will be allowed to depend on the 
latitude 
| p=p(9), | | 
because this greater generality willl be required in the physical discussion 
at the end of §3. By conditions (i)-(iv), the conservation equations 
become simply 
E?2- sin? 0(H°3+ pw?y?)=0 
tan 6 0H??/d6+ H?2—sin? 0(£33+- pw?y?)=0 

tan 0 0H3/00+ 3E%—0 
tan 0 0H?4/00+ E74=0 


since it is evident from the conditions of the problem that H*”’ is 


independent of ¢ and ¢. 
Addition of the first two of eqns. (12) leads to 


ea a2) 


a (sin? 6B22)=0, 


whence it follows from condition (iv) and the remark about it just after 
eqn. (11) above, that 


B2=0. MSTA! oc Seca wan TS) 
From (13) and the first of (12), 
E3= — pw*y?, ean WR ewe Ee) 
and by (6) 
peter nea? winltGy?;0 ts). (15) 
E44=pw'a‘ sin‘ 6y?. earner mas bes (16) 


The third and last of (12) are equivalent to 


0 jan8 A F128) — a og 24)__¢) 
ap (sin OE )=0, ap (sin bE ) 2 


854. L. Bass and F. A. E. Pirani on the 


which, by an argument analogous to that resulting in (13), lead to 
E®=0, oe Ee et 


Ha: wha) 20g inn & 
Equation (5) implies that 
URE he N HED sia. Th sil ls OD 
Thus the elastic tensor #” is prescribed completely. Rejecting all terms 
of order higher than the second in w, one is left with only one non- 
vanishing component of H“’, namely H*%, given by (14). The only 
non-vanishing components of the total energy-momentum tensor 7” are 
now found from (3) to be 
fh tom ive rip We Se Spy eo See hipPo noel See) 


Transformation to rectangular coordinates 


x=r sin 6 cos ¢, y=r sin 6 sin ¢, z=r cos 6 
yields 
0, 0, 0, twa sin # sin ¢ 
0, 0, 0, —twa sin 6 cos d 
T= OX: . 
J 0, 0, 0, a) 
twa sin @ sing, —twasin@cos¢, 0, 1 


(21) 
Comparison with Thirring’s expression for 7',,, (1918, eqn. (6)) shows that: 
in the present reference frame, all components of the total stress vanish 
upon the introduction of H”. The elastic stresses are just balanced by 
the centrifugal forces acting on the shell. 


§ 3. THe APPROXIMATE INTEGRATION OF THE FreLD EQUATIONS 


In this section the influence of the rotating shell on the metric near 
its centre will be determined, up to order w? in the angular velocity. 
We shall only outline the calculations, which follow closely those in 
Thirring’s papers (1918, 1921), where the details may be found. Let the 
linear approximation to the field be given by 


G9= —S yw TEV py | 
where (22) 
VV a BOs Vine | 
uy py 2° av / aM 
and 
, k T (ZL, Y,%; t—R) ‘ 
Vy =— 16,3 Seah Lia ls R dV 4. ni es Wee (ad 


Here & is the Newtonian constant of gravitation and R the distance 
from the (3-dimensional) volume element: dV,=a* da sin 6 dé dé to the 
point P near the centre of the shell where Yu 18 to be evaluated. The 
coordinates of P shall be 


*=rsin 6, y=0, z2=r cos 6. 
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The expression R~1 will be expanded in powers of ra~1, powers higher 
than the second being neglected (cf. Thirring 1918, eqn. (9)). Thus 


Rk-1=a-!{1+ra-\(sin 6 cos ¢ sin 6)-+-cos 6 cos 6) —4r2a-2 
+ $r?a-*(sin 6 cos ¢ sin 6)+cos 6 cos 6,)?}. Wes... (24) 


Besides adding the elastic term to 7”, we shall generalize Thirring’s 
work in the following respect: we shall assume that the density of 
the material of the shell varies with the latitude according to 

p=p)(1+Nw?a? sin? 6), lise ee het heel 25 ) 
where py and N are constants. Our purpose is to investigate the effects 
that would result from compensating the special-relativistic increase in 
mass of the moving shell. By a straightforward calculation of (23), 
from (21), (24) and (25), one finds for the components of the metric tensor 
near the centre of the shell (after a rotation about the z-axis to a general 
point P) 


mad vag 0 0 sik Mwa- ly 
il aya 0 Sik Moat 
Ke 0 0 eres yet 0 : 
SikMaa'y = —4ikMwa-'x 0 —1+37a0 
(26) 
where 
Yaa = 4k Ma-1{1 + $0202 + E(1+ N) w(x? + y?— 22) }, 
while 


M=| p dV y=4npy? da{l+4(1+2N)a%?} Lnien 


is the total mass of the shell. 
The equations of a geodesic near the centre of the shell reduce for 


small velocities to 
t= —8&kMoway+2(1+N)kMwra- tx | 
j= *kMoaez+2(1+N)kMora-ty 7 ee pr) 
— —s(1+N)kMo"a~tz j 


The ‘Coriolis’ force (w-terms) is naturally the same as that found by 
Thirring, since all the additional terms which we have introduced are 
of order w?. The ‘centrifugal’ force (#?-terms) has, compared with 
Thirring’s, an additional factor 4(1+N). Thus for N=0 (the case 
considered by Thirring), the centrifugal force is just halved by the 
introduction of the stresses in the shell. The choice N=—3, which, 
as may be seen from (27) corresponds to a shell with total mass equal 
to that of a shell at rest with uniform density po, reduces the centrifugal 
force to one quarter of Thirring’s result. 
The most interesting case is 


N=—1. 
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This corresponds to a mass distribution which just compensates for the 
special-relativistic increase in density, and represents a uniform mass 
distribution in the reference frame in which the shell is rotating. This 
choice leads, by (28), to the annihilation of both the radial and the axial 
‘centrifugal’ forces, leaving the ‘ Coriolis’ force intact. 

Thus if Thirring’s assumptions are modified merely by the introduction 
of the elastic stress term (which is essential if the conservation law is 
to be satisfied) the surviving w?-terms in (28) may be attributed entirely 
to the variation in density between the poles of the shell and its equator 
(a consequence of the special-relativistic mass increase with motion) and 
may be abolished by giving the rotating sphere a uniform mass distribution. 
This makes the signs of the w?-terms immediately plausible, as well as 
their independence of the sense of rotation. 

If this clears up the appearance of an axial component (towards the 
centre) of the centrifugal force, unexpected and unwanted from the point 
of view of Mach’s principle, it throws further doubt on the interpretation 
of the w?-terms as an ‘induced centrifugal force’ as demanded by this 
principle (cf. Einstein 1945). We do not propose any resolution of this 
problem. If Newton’s bucket of water rotates, and the stars are at rest, 
then the water rises up the sides of the bucket to a paraboloid. According 
to Mach, this situation is indistinguishable from that at which the bucket 
is at rest and the stars rotate about it. But in the case of the rotating 
shell problem the boundary conditions are always Minkowskian—only 
the shell spins, and not the ‘distant’ stars. 
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ABSTRACT 


An interferometer, operating on the rotating-lobe principle, is described. 
It is shown how the frequency of the interference patterns observed from 
discrete radio sources may be controlled and made independent of the 
baseline length and direction. The application of the instrument to 
observations with a short north-south baseline and a long (630 A) baseline 
is illustrated by experimental results. It is pointed out that the use of 
this technique is of particular value for, the observation of weak sources 
with long baselines. 


§ 1. LyrRopUCTION 


AN important experimental problem in radio astronomy is to measure the 
apparent angular diameter and the distribution of intensity across the 
discrete sources. In general this measurement demands a high resolving 
power which can only be obtained by the use of an interferometer. It is 
the purpose of the present paper to describe an instrument which has been 
developed for this type of work and which has been termed a rotating- 
lobe interferometer. It is shown that the use of this instrument offers 
certain advantages over previous methods. 


§2. Tur PrincrPLe oF A RotTatine-LOBE INTERFEROMETER 


The principle of a rotating-lobe interferometer can be explained simply 
in terms of the effective polar diagram of the aerial system. Figure 1 
represents two aerials A and B which are separated by a baseline of length 
d and are connected by identical cables to a multiplier. In one of the two 
cables there is a continuously variable phase-shifter. If radiation from a 
distant point source is incident on the aerial system at an angle 4, to the 
normal to the baseline, for any given value of the phase-shift 4, the time- 
average of the multiplier output will be proportional to 


2rd 


P,, [9.4(9) 9p(@o) 12” cos | 1 sin 6y+4 | Pemaarce ss (1) 


where P, represents the flux from the source and g 4(99), gg(4) represent 
the power gains of the two aerials. The response of the whole system to a 


* Communicated by Professor A. C. B. Lovell. 


858 R. Hanbury Brown et al. on a Rotating-Lobe 


point source at any angle may therefore be represented as a polar diagram 
of lobes of alternate polarity as shown in fig. 2, the angular position of the 
central lobe of the pattern with respect to the normal being determined by 
the value of the phase-shift ¢. 


Fig. 1 


MULTIPLIER 


RECORDER 


A simple form of interferometer. 


Fig. 2 


a Sek Bc ln aA A Sa 
The effective polar diagram of reception of the simple interferometer shown in 
fig. 15 


As a discrete source transits the normal to the baseline, it will produce an 
interference pattern on the recorder with an apparent frequency f, given 


by 1 
fey, | on ; 008 3.608 «. | Sere not aay iL: 
“7 


where wz radians sec~ is the rate of angular rotation of the earth, A is the 
wavelength, 5 is the declination of the source and « is the angle between 
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the direction of the base-line and a line east-west. If the value of ¢ is now 
made to vary slowly with time at a rate w, radians sec~!, then the lobes 
of the polar diagram will appear to rotate, while the envelope of their 
maxima remains stationary. Under these conditions the frequency of the 
_ output pattern is modified and becomes 


tae d 
I= 5 ops 0085 008 1-60 | een ore ey (3) 


where the choice of sign depends upon the direction of rotation of the 
phase-shifter. It is therefore possible, by controlling both the speed and 
direction of the phase-shifter, to vary at will the frequency of the output 
pattern. Thus, by rotating the lobes in the same direction as the apparent 
angular motion of the source, the frequency of the pattern may be 
reduced indefinitely, while by rotating the lobes in the opposite direction, 
or more rapidly in the same direction, the frequency can be increased to 
any desired value. As discussed later this facility is of particular value 
when it is necessary to operate an interferometer with very short or very 
long baselines, or when the baselines are in a direction close to north— 
south. 


§ 3. DESCRIPTION OF A PRACTICAL EQUIPMENT 


3.1. A Simplified Outline 


A rotating-lobe interferometer can be. constructed from a number of 
different circuit arrangements, since there are many ways* in which 
the relative phase of the signals from the two aerials may be varied. The 
choice of any particular arrangement must necessarily depend upon the 
use for which the instrument is intended. The instrument described here 
has been designed specifically to measure the distribution of intensity 
across both weak and intense sources, using baselines which may be of 
widely different lengths and oriented in any direction. 

A simplified block diagram of the instrument is shown in fig. 3. A and 
B represent the two aerials which are connected by cables to separate 
receiving channels. The output from aerial A is amplified at the signal- 
frequency f, and fed to a mixer, where it is heterodyned by a local oscil- 
lator of frequency f,. The output from aerial B is treated in the same way, 
except that the frequency of the local oscillator is f, and therefore differs 
from that used in the other channel. The outputs of the two mixers are 
amplified at intermediate-frequencies (f)—/,), (fo—/2), and are combined 
in a detector to form a beat-frequency (f,—f.). This beat-frequency is 
amplified and rectified in a phase-sensitive rectifier. ‘The reference voltage 
of the phase-sensitive rectifier is generated by combining the two local 


* An interferometer using a rotating phase-shifter in one aerial, but no phase- 
sensitive rectifier, has been described by Little and Payne-Scott (1951). Also 
Bolton and Slee (1954) have used an interferometer in which the relative phase 
of the local oscillator at two mixers was varied slowly to produce a moving 
lobe-pattern. 
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oscillators in a mixer ; the beat-frequency between the oscillators is then 
fed through a continuously variable phase-shifter to operate the rectifier. 
The output of the phase-sensitive rectifier is passed through a time- 
constant and is displayed in the usual way on a pen recorder. 

Following the method used above it is convenient to discuss the operation 
of this system in terms of its effective polar diagram of reception. The use 
of two local oscillators of different frequency may be regarded as a‘method 
of introducing a continuous variation of 27 (f,—/,) radians sec~+ into the 
relative phase of the signals from the two aerials. Thus the lobes of the 


Fig. 3 


ne 
AMPLIFIER. 


A simplified block diagram of the interferometer. 


polar diagram appear to rotate at A/d (f,—/,) radians sec-!, and a discrete 
source in transit will produce an alternating component in the output of the 
detector given by 


P,[94(%) 9(0o) G4 Gp]? cos E A . Cos 6 Cos aE 2a (fo) | sat (a) 


where G4, Gz are the power gains of the two channels. 

The output of the phase-sensitive rectifier is proportional to the ampli- 
tude of the input wave and to the cosine of the phase-angle between this 
wave and the reference voltage. It therefore consists of an alternating 
component of low-frequency given by 


d 
PLa4(0) 930) Ga Ga)? 008| ops cos dcosztuy | . . (6) 


where w, is the rate of change of the phase of the reference voltage and 
the choice of sign depends upon the direction of the phase-shifter. 
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A comparison of eqns. (4) and (5) shows that the use of a phase-sensitive: 
rectifier brings the rotating lobes apparently to rest, and that the fre- 
quency of the interference pattern is determined in the usual way by the 
rotation of the earth if the phase-shifter is stationary. If, however, the 
_ phase-shifter is varied, then the angular frequency of the pattern is changed 
by the same amount. Thus it follows from eqns (3) and (5) that the 
phase-shifter shown in fig. 3 may be used to control the frequency of the 
interference pattern in precisely the same way as a phase-shifter is one 
aerial cable (fig. 1). 


3.2. Some Details of the Apparatus 


A schematic diagram of the complete apparatus is shown in fig. 4; 
although it differs in detail from the simple scheme in fig. 3, the principle 


Fig. 4 


"FIXED STATION MOBILE STATION 


A complete block diagram of the interferometer. __ 

(1) Preamplifier (158-2 Mc/s) ; (2) First mixer ; (3) Frequency multiplier { x 12) 
(4) Intermediate frequency amplifier (10 Mc/s) ; (5) Amplifier (14-02 Mc/s) ; 
(6) Automatic gain control; (7) Second mixer; (8) Intermediate fre- 
quency amplifier (10 Mc/s); (9) Intermediate frequency amplifier 
(3°76 Mc/s); (10) Delay cable; (11) Crystal oscillator (6-240 Mc/s) ; 
(12) Crystal oscillator (6-241 Mc/s); (13) Second mixer; (14) Auto- 
matic gain control; (15) Intermediate frequency amplifier (3-76 Mc/s) ; 
(16) Mixer; (17) Intermediate frequency amplifier (3-76 Me/s) ; (18) 

- Magslip phase-shifter driven by velodyne ; (19) Multiplier ; (20) Selective 
amplifier (1000 c/s) ; (21) Phase-sensitive rectifier ; (22) Time-constant ; 
(23) Recorder ; (24) Preamplifier (158-2 Me/s) ; (25) First mixer ; (26) 
Frequency multiplier (x12); (27) Crystal oscillator (14-02 Mc/s) ; 
(28) Intermediate frequency amplifier (10 Mc/s); (29) Amplifier 
(14-02 Mc/s). 
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of operation is unchanged. ‘The equipment is divided into a fixed and a 
mobile station. The signal received at the mobile station is heterodyned 
to a convenient intermediate frequency (10 Mc/s) and is transmitted over a 
coaxial cable together with the fundamental frequency of the local oscil- 
lator. At the fixed station this signal is amplified at 10 Mc/s and fed to a 
second mixer ; the local-oscillator fundamental is multiplied and used to 
heterodyne the signal received by the aerial at the fixed station. This 
latter signal is also amplified at 10 Mc/s and fed to a second mixer. The 
frequency difference between the two channels is introduced at the second 
mixers by the use of two second local oscillators which differ in frequency 
by about 1000 c/s. The signals are then amplified at the second inter- 
mediate frequency and are combined in a multiplier. The signal from the 
fixed aerial is delayed in the second intermediate frequency amplifier to 
compensate for the delay in the signal from the mobile aerial. The output 
of the multiplier contains a beat-frequency of 1000 c/s which is amplified 
in a selective amplifier and applied to a phase-sensitive rectifier. The 
beat-frequency between the two second local-oscillators is formed in a 
separate mixer and is fed to the phase-sensitive rectifier through a 
rotating magslip phase-shifter driven at an adjustable speed by a velodyne 
motor. 


§ 4. THE EQUIPMENT IN OPERATION 


The instrument described above was designed to measure the apparent 
angular diameter of the discrete sources observed in a previous survey 
with the large fixed paraboloid (218 ft.) at Jodrell Bank. For the present 
purpose the large paraboloid was used as the fixed aerial of the interfero- 
meter and the mobile aerial usually consisted of a small array of about 
35 m? aperture ; however, for observations of the most powerful sources 
the mobile aerial was occasionally a single dipole. The measurements 
were carried out by observing the sources with baselines of different length 
and direction as they transited the meridian, and their apparent angular 
diameter was found from the well-known relation between the distribution 
of intensity across a source and the relative amplitude of the interference 
pattern at different aerial spacings. 

The results of the first series of measurements have already been 
reported (Hanbury Brown, Palmer and Thompson 1954). It was found 
that several of the intense sources have angular diameters of the order of 
one or two degrees, hence these sources were completely resolved with 
baselines of less than about 60. On the other hand a number of the 
sources, particularly those at high galactic latitudes, showed little trace 
of being resolved with a baseline of 500 A. 

These preliminary measurements illustrate clearly the need for an 
interferometer which will operate over a wide variety of baseline lengths ; 
furthermore since the sources may not be radially symmetrical, as in the 
case of the intense source in Cygnus, it is essential that the instrument 
‘should operate with a baseline in any direction. Equation (2) shows that, 
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if the frequency of the pattern is determined only by the rotation of the 
earth, then it is proportional to the length of the baseline and to the cosine 
of the angle between the baseline and a line east-west. Thus for very short 
baselines, or for baselines close to north-south, the frequency of the 
pattern is low. If the baseline is comparable in length with the aperture 
of either. of the aerials, it becomes difficult in practice to measure the 
amplitude of the interference pattern during the transit of the source 


Fig. 5 


—20MINs. —10 mins. TRANSIT +10 MINs. + 20 Mins. + 30mins. 


5 (a) 


5(b) 


The use of a rotating lobe interferometer to decrease the frequency of an 
interference pattern. 

(a) Record of a weak source observed with a baseline of 630 A in a direc- 
tion east and west, taken with the phase-shifter stationary. The 
output time constant is 4 seconds. 

(b) Record of the same source taken under identical conditions but with 

the phase-shifter rotating and decreasing the frequency of the pattern, 
and the output time constant increased to 30 seconds. 
N.B. For both these records the 218 ft. paraboloid was used as one 
aerial of the interferometer. At a frequency of 158-2 Me/s 
the beam is narrow and the interference pattern is therefore 
restricted to about 10 minutes on either side of transit. 


through the aerial beam. Conversely, if the baseline is long, the frequency 
of the pattern becomes inconveniently high and it is necessary to reduce 
the time-constant of the receiver output. Since the signal-to-noise is 
proportional to the square-root of the time-constant, it follows that, for 
long baselines, the overall sensitivity of the equipment will be reduced. 
Practical experience shows that the difficulties outlined above can be 
substantially reduced by the use of the equipment described in the present 
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paper. Observations* have been made with baselines in several directions 
and with a variety of lengths between 21 \ and 630 A. It was found that, 
with the exception of very short baselines as mentioned later, it was 
always possible to bring the frequency of the interference pattern to a 
convenient standard value by adjusting the speed and direction of the 
phase-shifter. As an example fig. 5 shows a record of a weak source 


Fig. 6 
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The use of a rotating lobe interferometer to increase the frequency of an inter- 
ference pattern. 


(a) A record of an intense source taken with a baseline 50 A in a direction 
north-south and with the phase-shifter at rest in an arbitrary position. 


(6) A record taken under identical conditions but with the phase-shifter 
rotating and increasing the frequency of the interference pattern. 


observed with a baseline of 480 A and with the phase-shifter stopped ; 
while fig. 5 (b) shows a record of the same source, taken under identical 
conditions, but with the phase-shifter adjusted to slow down the inter- 
ference pattern to a convenient value, and a greater output time constant. 
It is clear that a considerable increase in the effective signal-to-noise ratio 
has been gained by the use of this technique. 

As an example of the reverse operation, namely the speeding up of a 
pattern, fig. 6 (a) shows a record of the transit of the intense source in 
Cassiopeia made with a baseline of 50 A oriented north-south and with the 
phase-shifter at rest in an arbitrary position. The amplitude of this 
pattern varies with the position of the phase-shifter and in the example 
shown it happens to be small. Figure 6 (b) shows a record of the same 
source made under identical conditions, but with the phase-shifter rotating 
at a speed chosen to give a pattern of convenient frequency. 


* For observations with baselines exceeding 500 A the cable was replaced by 
a radio link. 
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Although the equipment operated satisfactorily under nearly all 
conditions, some difficulty was experienced when the phase-shifter was 
used to increase the frequency of patterns from weak sources observed 
with very short baselines. Thus, when the aerials were separated by a 
distance of less than about 40 A, approximately the diameter of the large 
paraboloid, a small and almost constant output was observed from the 
multiplier even in the absence of a source. If the phase-shifter was 
stationary, this output represented a small displacement of the zero of the 
records and was unimportant ; however, with the phase-shifter rotating, 
the output was converted to a sine-wave and tended to confuse the 
records of weak sources. It is to be expected that, when the baseline is 
comparable with the dimensions of the aerials, there should be an output 
due partly to receiver noise cross-coupled between the aerials and partly 
to the background radiation from the sky. This problem has not yet been 
investigated completely ; however it appears that in the present equip- 
ment the limitation is set by the considerable mutual impedance which 
exists, at close aerial spacings, between the mobile aerial and the elevated 
primary feed of the fixed paraboloid. The extent to which the perfor- 
mance at short baselines can be improved by reduction of this mutual 
impedance remains to be investigated. 


§ 5. Discussion 

In the simplest form of interferometer the output powers of two 
spaced aerials are added together, either before or after amplification, 
and the total power is recorded. Thus, in addition to the interference 
pattern produced by a discrete source, the output contains components 
proportional to the receiver noise and to the background radiation. The 
presence of these large components imposes serious limitations on the 
observation of weak sources and a considerable gain in effective sensitivity 
results from their removal. In principle this can be accomplished by the 
use of the simple system shown in fig. 1 in which the outputs of the two 
aerials are multiplied rather than added. But in practice this method 
is not attractive since the design of a satisfactory multiplier is difficult 
when the signal is very small compared with the noise. 

An effective way of overcoming this difficulty has been developed by 
Ryle (1952) and is termed ‘ phase-switching ’. In this system the phase 
of one of the spaced aerials is reversed periodically by a mechanical 
switch and therefore, when the outputs of the two aerials are combined, 
signals which are mutually coherent in the two aerials produce switch- 
frequency components in the amplitude of the combined output. These 
components are selected from the output of the receiver and are rectified 
in a phase-sensitive rectifier operated in synchronism with the phase- 
reversing switch. It has been shown (Ryle 1952) that this system removes 
the mean value of the receiver noise and background radiation from the 
output of the phase-sensitive rectifier. It can also be shown that the 
effective polar diagram of the aerial system is a series of alternate positive 
and negative lobes similar to that shown in fig. 2. 
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The interferometer described in the present paper also eliminates the 
mean value of the unwanted components from the output. In fact the 
principles of operation of the present system and the phase-switching system 
are similar, and the latter may be regarded as a special form of rotating- 
lobe interferometer in which the lobes are rotated in discrete jumps rather 
than continuously. A comparison of the two systems suggests that, in 
general, the phase-switching technique is preferable for the precise 
measurement of the celestial coordinates of discrete sources. For this 
type of measurement it is essential to know the relative phase of the 
signals in the two aerials and, if the baseline is not too long, this can be 
best achieved by connecting the aerials to a common receiver through 
identical cables. Under these conditions phase-switching has the con- 
siderable practical advantage that it is much simpler to introduce a 
phase-reversing switch into one of the cables than a continuously variable 
phase-shifter. On the other hand, if the baseline is so long that separate 
amplifiers must be used for each aerial, or if a precise knowledge of the 
relative phase of the signals is unimportant, then the technique described 
in the present paper has certain advantages. For example the use of two 
local oscillators of different frequency is reasonably simple and avoids the 
use of a mechanical phase-reversing switch ; furthermore the frequency 
difference between these oscillators, which corresponds to the phase- 
switching frequency, is not limited by mechanical considerations. 

The most interesting property of the present instrument, and one which 
offers the greatest potential advantage over previous methods, is the 
simplicity with which the frequency of the interference patterns can be 
controlled. This property is likely to prove useful for interferometers 
working on long wavelengths with restricted baselines, for the measurement 
of the intensity distribution across sources using baselines close to north— 
south, and for observations of weak sources, or the detailed structure of 
intense sources, with long baselines. The last of these applications, 
namely the use of long baselines, is of particular interest. Thus it follows 
from the analysis given above that, by the proper use of lobe-rotation 
technique, the resolving power of an interferometer can be increased 
without loss of sensitivity until the limit is set by other factors such as the 
instability of propagation. 
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ABSTRACT 


The surface energies of the alkali metals sodium, potassium and 
lithium and their temperature coefficients have been measured experi- 
mentally. The values for lithium and sodium are in good agreement 
with similar theoretically calculated surface energies obtained by Stratton 
(1953). The disparity between the experimental and theoretical surface 
energy values for potassium and that between the experimental and 
theoretical temperature coefficients of the three metals is discussed. - 


§ 1. LIyrRopuUcTION 


A NUMBER of theoretical approaches have been made to the calculation 
of the surface energies of metals, but the final verification of such treat- 
ments was not possible due to the lack of accurately measured experimental 
values. In particular, experimentally determined surface energies for the 
metals lithium, sodium and potassium are germane since the physical 
basis of the majority of these calculations is: most nearly fulfilled by 
these metals. 

Prior to this study, experimental data on the surface energies of the 
three metals was incomplete. No measurements had been made on 
lithium and such values as existed for potassium (Quarterman and 
Primak 1950) were neither accurate nor comprehensive... The only 
accurate data for sodium was that reported by Poindexter and Kernaghan 
(1929) which leads to a surface energy value of 206 ergs/cm? at 98°c and 
a temperature coefficient of —0-05 ergs/em?/°c. A second determination 
on sodium has recently been reported by Addison, Kerridge and Lewis 
(1954) ; their results are referred to subsequently. 

The surface energies of the metals sodium, potassium and lithium 
have now been determined experimentally and these values allow a 
more exact comparison with calculated values, thus leading to a verifi- 
cation of the analytical methods employed. 


§ 2. EXPERIMENTAL 


The maximum bubble pressure method was employed in the experi- 
mental study for the following reasons. A fresh metal surface is formed 
in successive determinations thus minimizing surface-contamination 
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errors, always a danger with the alkali metals. The measurements may 
be effected remotely with no direct observation on the liquid metal being 
necessary. The physical parameters involved are well defined, while an 
analysis by Porter (1930) has established the mechanics of bubble 
formation and detachment. Finally there is a direct theoretical relation- 
ship between the parameters measured and the surface energy of the 
liquid; this obviates any need for an empirical calibration procedure. 

The maximum pressure required to form and detach a bubble from 
an orifice immersed in a liquid may be analysed as follows : 


20 
Prax=9-h.(p— po) = ; “he, eh Mle ROD 


where P,,,,—maximum pressure, h=depth of immersion of the orifice tip, 
p=liquid density, pp=gas-phase density, c—surface energy of the liquid, 
and r=effective orifice radius. 

In applying eqn. (1) it is necessary to ascertain whether 7 is the 
external, internal, or some intermediate radius of the orifice at the point 
of detachment of the bubble. Furthermore, a correction must be applied 
for the distortion of the bubble from a spherical form due to the gravita- 
tional effect of the varying head of liquid within it. Porter’s analysis 
(1930) was employed in the selection of the effective orifice radius while 
Sugden’s method of correction (1922) was applied for bubble distortion ; 
both factors have been discussed at length by the author (Taylor 1954). 

The metals employed were prepared by vacuum distillation under 
conditions designed to remove both heavy metal and metalloid impurities ; 
analyses of the metals are contained in table 1. 


Table 1. Analyses of Alkali Metals 


Metal Purity (%) | Impurity content (p. p. m.) 


Na 99-995 B, 40; K, 20; O, 5-10 
K 99-895 Na, 10 
Li >99-95 Na, <100; Cl, <400 


Steps were taken to ensure that this purity was maintained during the 
transfer of the metals from the purification still to the measuring 
apparatus. 

This apparatus is shown diagrammatically in fig. 1. The orifice tube (C) 
is movable, via a Wilson seal, through the cover plate of vessel (A), the 
positioning of the orifice tube being facilitated by means of the circuit (H). 
The latter is ‘made’ when the orifice tip first touches the liquid metal 
surface and the levelling arms (G) and (G,) may then be zeroed. Sub- 
sequently the orifice tip may be located at any required depth by means 
of a screw-gauge applied over the levelling arms; the accuracy of this 
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setting is +0-0005cm. Orifice tubes were made of mild steel for 
measurements in lithium and of stainless steel for sodium and potassium. 
The orifice tip was given a metallographic polish, on a surface perpendicular 
to the tube axis; the tube diameters were measured optically to an 
accuracy of +0-01%. The manometer (E) could be read under ideal 
conditions to an accuracy of +-0-02 cm but in practice a slight pulsing” 
of the manometer levels caused this error to be slightly higher; low 
density, low vapour pressure vacuum oil was employed as manometer 


Fig. 1 
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Apparatus (diagrammatic). 


fluid. Bubble formation and detachment is effected through the tombac 
bellows (F), a fine screw-thread attached to the latter ensuring a slow 
and controlled rate of bubble formation. Cleanliness of the apparatus 
was ensured by prior electropolishing of all metal parts coming in direct 
contact with the liquid metals, while a high-vacuum system ensured 
freedom from gaseous contamination ; high-purity argon was employed 
as the blanket gas during the measurements. 

The transfer of the high-purity metals into the apparatus was effected 
via a steel bulb which was filled directly from the purification still, and 
attached subsequently to the measuring apparatus. No contamination 
of the liquid occurred during this transfer process. 
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Surface energy measurements were normally made on a charge of 
150-200 cubic centimetres of liquid at three separate depths at any one 
temperature. The average of six separate determinations was taken as 
the maximum pressure at any given depth and temperature. 


§ 3. RESULTS 


The surface energy—temperature relationships determined experimentally 
are presented in figs. 2, 3 and 4 for sodium, potassium and lithium respec- 
tively ; the experimental data were corrected for bubble distortion in 
accordance with Sugden’s method (1922). The surface energy values at 
the respective melting temperatures are: 191 ergs/em? for sodium (melting 
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Surface energy—-temperature relationship for sodium. 


temperature 98°C), 101 ergs/cm? for potassium (melting temperature 64°0) 
and 398 ergs/cm? for lithium (melting temperature 179°C), while the corre- 
sponding temperature coefficients of surface energy are :—0-1 ergs/em?/°o, 
—0-11 ergs/em?/°c and —0-14 ergs/em?2/°c. 
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Fig. 3 
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Fig. 4 
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Surface energy—-temperature relationship for lithium. 


The accuracy of these experimental surface energy values is controlled 
by the errors involved in the several physical parameters employed in. 
the final calculations. These sources of error are listed in table 2. 
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Table 2. Percentage Error in Physical Parameters 


Metal Na K Li 
Physical parameter Percentage error 
Density of manometer fluid +0:1 +0:1 +0°1 
Orifice diameter +0-01 +0-01 +0-01 
Orifice depth setting 7 +0-1 +0-1 +0-1 
Density of liquid metal +0-1 +0-1 +0-2 
Maximum pressure measurement +2-0 4-2-5 +0°3 


In the case of sodium and potassium the maximum error is in the 
individual maximum-pressure readings on account of the low surface 
energy component of the maximum pressure for these liquids; this 
error is relatively smaller for lithium due to its higher surface energy. 
The fact that pressure readings were made at three separate depths at 
any one temperature probably reduces the magnitude of this error 
-slightly and it is considered that the values are accurate to +2° in the 
case of sodium and potassium and +1°, for lithium. This is the order 
of magnitude of the scatter of the experimental points shown in figs. 2, 3 
and 4. In addition, density—temperature relationships, derived from the 
measurements at different depths, show a similar accuracy in relation to 
previously determined values (fig. 5). 
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Density-temperature relationships for sodium, potassium and lithium. 
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Finally the present experimental data for sodium is in good agreement 
with a recent independent determination (Addison, Kerridge and Lewis 
1954) which gives a surface energy of 195 ergs/em? at 98°c and a 
_ temperature coefficient of —0-1 ergs/em2/°c. 


§ 4. THEORETICAL CALCULATIONS OF SURFACE ENERGY 


A number of theoretical calculations of the surface energies of metals 
have been made and these have been discussed in detail by Ewald and 
Juretschke (1953). Hence only a brief reference will be made to the 
several treatments with which the present experimental data are compared. 

Two of the earliest calculations were those by Gogate and Kothari (1935) 
and by Samoilovich (1945). Employing a free-electron model the former 
considered the electron motion in the surface of a metal which, for this 
purpose, was regarded as a two-dimensional plane film. From a con- 
sideration of the kinetic energy of electrons confined within such a surface 
film, a relationship was established between the total energy, the kinetic 
energy, and the surface energy of the system in relation to its surface area. 
The treatment by Samoilovich (1945) was based on the density gradient 
which exists at a metal surface. An anisotropic strain tensor is thus 
produced at the surface so that, when the surface area is increased, energy 
is expended in doing so. The overall force-density within the metal was 
calculated using the Thomas—Fermi—Weizsicker approximation (1935) 
and this density was then equated to that in the surface. 

Fundamentally the two foregoing analyses cannot be regarded as 
accurate due to the assumptions made regarding the boundary conditions 
at the surface of a metal. Furthermore the accuracy of the Weizsacker 
approximation appears to be in some doubt (Ewald and Juretschke 1953) 
and leads to excessively high surface energy values. 

Subsequent treatments have given more detailed consideration to the 
exact boundary conditions existing at a metal surface. In particular, 
the existence of a definite potential step at the surface is recognized 
with its consequent effect on electron distribution. Thus Brager and 
Schuchowitzky (1946) calculated, by a wave-mechanical treatment, the 
energy increase associated with the formation of a new surface, assuming 
an infinite potential step at the surface. However, this step at the 
surface of a metal has a finite value and for this reason the previous 
treatment is oversimplified. Huang and Wyllie (1949) have extended 
the analysis to the case where the potential step has a finite value. In 
addition they have included a term in the expression for surface energy 
which takes account of the electrostatic charge at the metal surface 
consequent on the irregular charge distribution in this region. 

Further refinements of the surface energy calculations have been 
proposed by Huntington (1951) and also by Stratton (1953), both of 
whom consider the phenomenon of charge conservation at the metal 
surface. Thus Huang and Wyllie (1949) defined both the position and 
the height of the potential step and consequently it was uncertain 
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whether the negative charge external to the surface would be completely 
balanced by the positive charge within. Of the two modified treatments, 
that by Stratton (1953) is the more complete. The analysis employed 
a free-electron gas model and established the surface energy component 
of the total free energy of the metal by considering the energy of the 
electron gas confined within a finite potential well. The electronic charge 
density remote from the surface was balanced by the positive ionic charge 
density inside the potential well. The final surface energy expression 
included a contribution from the electrostatic surface double-layer and 
a term correcting the surface energy at absolute zero for effects at higher 
temperatures. 


§ 5. COMPARISON OF THEORETICAL AND EXPERIMENTAL VALUES 


The several calculated surface energies of the metals sodium, potassium 
and lithium are compared with the experimental values in table 3; 
where the theoretical treatments permit, the values are compared at 
the respective melting temperatures. Temperature coefficients of the 
surface energies are also compared where possible. 


Table 3. Calculated and Experimental Surface Energy Values 


Surface energy 
at the melting 
temperature 


Temperature 
coefficient of 
surface energy 


Author(s) (ergs/em?) (ergs/em?/°c) 
Na K Li Na K Li 
Gogate and Kothari (1935) | 278 | 121 639 | —0-22 | —0-24| —0-18 
Samoilovich (1945) 400* | 224* | 655* | —* —* —* 
Brager and Schuchowitzk y 
(1946) 780% |) 300" Telio0<4) —— —* —* 


Huang and Wyllie (1949) | 440 L80 890 | —0-12 | —0-08 | —0-19 


Stratton (1953) 190 70 400 —0-06 | —0-04 | -—0-07 


Experimental 191 101 398 | —0-10| —0-11] —0-14 


* Values at 0°K ; no temperature effect considered. 


The surface energy values as calculated by Gogate and Kothari (1935), 
Samoilovich (1945), Brager and Schuchowitzky (1946), and Huang and 
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Wyllie (1949), are from one and a half to three times the experimental 
values in the case of the three metals. 

_ Stratton’s treatment (1953) appears to give excellent agreement between 
the calculated and the experimental values for sodium and lithium. 
However it must be recalled that the two sets of results are not strictly 
comparable since the theoretical values are for the solid at the melting 
temperature, no term being included for the change in surface energy 
on melting. An order of magnitude to this change on melting may be 
obtained from the relative interfacial tensions between the solid and 
liquid phases at the melting temperatures; these interfacial tensions 
may be calculated from supercooling data (Turnbull 1950). On this 
basis the interfacial tensions were 30 ergs/cm? for lithium and 20 ergs/em2 
for sodium. Assuming that the difference between the surface energy of 
the solid and the liquid phase at the melting temperature is of the same 
order of magnitude as the interfacial tension, the corrected theoretical 
surface energy values would be 370 ergs/cm? for lithium and 170 ergs/cm? 
for sodium, values still in good agreement with the corresponding 
experimental data. 

Stratton’s calculated surface energy for potassium is considerably 
lower than the experimental value, and any correction to the calculated 
data for the change in surface energy on melting would only serve to 
accentuate the existing discrepancy. A factor accounting for this dis- 
parity may be the failure of the free-electron—gas model when applied 
to potassium. In Stratton’s analysis (1953) the periodic potential of 
the ionic core is replaced by an equivalent uniform potential field and 
this approximation is only true when the potential spread round the 
positive ion core is small. A similar effect appears to be reflected in 
Stratton’s values for the other monovalent metals copper, silver and gold 
where the calculated surface energy values are all considerably lower 
than the experimental values. 

In comparing the theoretical and experimental temperature coefficients 
of surface energy it appears that the best agreement was obtained in 
the case of Huang and Wyllie’s (1949) analysis. This agreement would 
seem to be somewhat fortuitous since their treatment of the effect of 
temperature on surface energy was oversimplified, that of Stratton (1953) 
being physically more exact. Stratton’s values are somewhat lower 
than those obtained experimentally but this difference may not be 
entirely due to errors in the two sets of values. Once again, the theoretical 
values are strictly applicable to the solid state only and in extending these 
to the liquid no account is taken of changes in lattice properties on melting. 
The changes in lattice properties which do occur on melting, e.g. increased 
lattice vibrations, are of such a nature that a higher temperature coefficient 
of surface energy would be expected in the liquid than in the solid phase. 
Thus the comparison of theoretical and experimental temperature 
coefficients of surface energies cannot be made rigorously. 
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§ 6. CONCLUSIONS 


The surface energies of the alkali metals sodium, potassium and lithium 
have been determined experimentally as: 191 ergs/em*, 101 ergs/em? and 
398 ergs/em? at the respective melting temperatures. The temperature 
coefficients of the surface energy were found to be :—0-10 ergs/cm?/°c for 
sodium, —0-11 ergs/em2/°c for potassium and —0-14 ergs/em?/°c for 
lithium. 

The most rigorous theoretical calculation (Stratton 1953) of the surface 
energies of metals leads to values for lithium and sodium which are in 
good agreement with the experimental data for these metals. The 
agreement between the experimental and calculated data is poor in the 
case of potassium but it is suggested that this may be a consequence of 
the failure of the free-electron—gas model when applied to potassium. 
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SUMMARY 


Ilford G5 emulsions exposed to a 950 Mev proton beam from the 
Birmingham synchrotron, have been searched for proton—proton elastic 
collisions. Forty-six such events have been found in 950 metres of proton 
track, equivalent to a cross section for elastic scattering of 15-5-+-2-5 mb. 
The angular distribution of the scattered protons in the centre-of-mass 
system shows a strong forward peak in agreement with results obtained 
at Brookhaven by cloud chamber and counter workers. 


; § 1. INTRODUCTION 

Unt1t recently studies of proton—proton scattering were limited to energies 
below 660 Mev. Proton synchrotrons have now extended the available 
range up to 1 Bev at Birmingham in Great Britain and to several Bev 
at Berkeley and Brookhaven in the United States. In this paper we 
report the preliminary results of an experiment to study the elastic scat- 
tering of 950 Mev protons by the hydrogen nuclei normally present in 
nuclear research emulsions. This technique, initiated by Goldhaber 
(1952, 1953), is being used by Kao and Clark (1954) to study proton—proton 
scattering at 420 Mev, and it has been widely used in the study of pion— 
proton scattering. Proton—proton scattering experiments at 900 Mev 
and 1000 mev have been briefly reported by the Brookhaven cloud chamber 
and counter groups respectively (Morris, Garrison, Fowler, Fowler, Shutt, 
Thorndike and Whittemore 1955, Smith, McReynolds and Snow 1955). 


§ 2. MeTHoD 


Ilford G5 normal nuclear research emulsions 2 in. x 2 in, x 400 w in size, 
were exposed to a single pulse of the scattered-out proton beam of 
the Birmingham synchrotron (Moon, Riddiford and Symonds 1955). The 
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energy of this beam is 950*3} Mev. The beam intensity used was about 
2 104 protons per square cm. The emulsions were searched by following 
each track along its length and all events found, including scatterings, 
were recorded. The low beam intensity reduced the chance of an observer | 
accidentally changing from one track to another during the scanning. An 
elastic interaction of an incident proton with a hydrogen nucleus in the 
emulsion will produce an event having two secondary tracks, i.e. an event of 
type 2p. To identify these interactions the projected angles, dips and grain 
densities of the primary and secondary tracks of every 2p event were 
measured, except when : 

(1) either track was identified by inspection as an «-particle or heavier 
nucleus, or 

(2) an additional track of length less than 5, due to a recoil nucleus, 
was observed. 

In order to calculate the true angles in space (say, «, and «,), of the 
secondary particles relative to the direction of motion of the primary particle 
it is necessary to know the shrinkage factor of the emulsion at the time at 
which the measurements are made. The method by which this was achieved 
has been described elsewhere (Duke, Lock, March and Munir 1955). 

When a moving proton collides elastically with a proton at rest the 
conservation laws yield the following equation for the relation between 
%1, %, and the total energy of the system H, measured in units of the proton 
rest-mass : 

tan «, tan «.=2/EH. Ce a eee 


For the energy 95013) Mev, tan «, tan «, will be 0-6679°?, Such an 
event must also be coplanar. This may be expressed by the condition 
that 6=0 where ¢ is the angle which the primary particle makes with the 
plane containing the two secondary particles. 


Now sees COS 5p COS 6, COS dy 
sin 0 


x [tan 6, sin 6,—tan 6, sin 0,+-tan 6) sin (8.—6,)]_ . (2) 


where 6, 5, and 6, are the dip angles of the primary and of the secondary 

_particles, 6, and 6, are the projected angles of the secondaries to the 
direction of motion of the primary and @ is the true angle in space between 
the secondary tracks. For an elastic collision between free protons 
@=«,-+%», which at this energy must lie between 78° and 90°. Events in 
which #<5° and @ is between 70° and 95° were carefully remeasured ; 
the errors in ¢ and © were usually less than 1°. The wide limits used 
ensured that no elastic interaction was lost. 

For each of these events we have inserted the larger of the two measured 
values of « into eqn. (1), and obtained the value expected for the other. 
The deviation 4, of the measured value of the smaller « from that 
expected should be zero for an elastic collision. Ona plot of 4 against 74 bi 
true elastic collisions will be represented by points within experimental 
error of the origin. 
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§ 3. REsuLts 

The path length followed in plates of known shrinkage factor is at 
present 710 metres and the total number of 2p events found is 390. 
Figure 1 is a plot of ¢ against 4, for those 62 events with <5° and 
|4,|<5°. There is a definite grouping of the experimental points around 
the origin. Consideration of the error involved in the individual deter- 
minations of ¢ and 4, leads us to expect that an elastic p—p collision will 
be represented by a point within a semi-circle of radius 1-:5°. Elastic 
collisions of the incident protons with protons bound in a nucleus probably 
account for the majority of the experimental points lying outside this 
boundary. Their distribution indicates that 2 of the 35 events within the 
semi-circle are also caused by such collisions. Thus we obtain for the 


Fig. 1 


A ee Ae Ce se ye ee 


A,, in degrees 


A plot of the degree of coplanarity, ¢, against the deviation from angular 
correlation 4, for 62 events of type 2p, with 6<5° and | 4,|<5°. 


mean free path, A, for elastic free proton—proton collisions the value 
710/33=21-6-+4-0 metres. Assuming the free hydrogen content of the 
emulsion to be (3-13--0-06) x 10? atoms per em? (Waller 1955) the cross 
section for elastic scattering is 148-25 mb. Systematic errors due to 
the selection procedure used are expected to be less than the quoted 
statistical uncertainty. 
In the preliminary stages of the experiment some scanning was carried 
out on emulsions whose shrinkage factor was not accurately known. Thus 
the values of 4, for the events found are liable to contain a rather large 
uncertainty, but the values of ¢ are almost independent of the shrinkage 
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factor. Following the same procedure as that outlined above, 13 elastic 
free p—p events were found in 237 metres of track, corresponding to a cross 
section of 17-+5mb. Forall 46 events in 947 metres A=20-6+3-5 metres 
and oc=15-5+2-5 mb. 

We may calculate the angles, «*, of the secondary particles with respect 
to the primary in the centre-of-mass (C) system by means of the relation- 
ship tan «=(1—f,2)/? tan «*/2 where f, is the velocity of the C-system. 
For an incident proton energy of 950 Mev f, is 0-58 and we have tan 
a*/2—1-23 tan «. The angular distribution thus obtained for the 35 
events within the 14° semi-circle together with the 13 events on the plates 
of unknown shrinkage factor is given in the table in units of 10° in orderthat 
it may be compared directly with the Brookhaven data, and in four 
equal units of solid angle. 


Anglein Number Morris etal. Fourequal Number  do/dw in mb 


C-system of normalized units of of per 
in degrees events - solid angle events steradian 
0 (24)-10 2 1 Q —414° 25 5-3+1-0 

10-20 8 9 414-60° 15 3-1+1-0 
20-30 7 ll 60 —754° 6 1-3+0°5 
30-40 6 ih 753-90° 2 0-4+0:3 
40-50 12 5 
50-60 5 4 
60-70 1 5 
70-80 6 3 
80-90 1 3 


{ The total number of events given by Morris et al. is 70. 


§ 4. Discussion 

The cross section for elastic p—p scattering which we have found is 
15-5+2-5 mb. A value of 21 mb may be obtained from the cloud chamber 
work of Morris e¢ al. using their value of 38-+-6 mb for the total cross section. 
If this latter quantity is taken to be 49-+-2 mb, as given by Shapiro, 
Leavitt and Chen (1954), then the data of Morris et al. yields an elastic 
cross section of 26 mb. Smith e¢ al. find a value of 19+3 mb, which is 
based on an extrapolation of their data from 30° to 0°. It appears 
possible that the elastic scattering cross section decreases slightly between 
400 and 1000 mev, for its value at 437 Mev is 23-8-+1-2mb (Sutton, 
Fields, Fox, Kane, Mott and Stallwood 1955). 

In the present experiment the minimum scattering angle which we can 
detect is approximately 1°, which corresponds to an angle in the O-system 
of 25°. The range of the recoil proton in sucha case is about 5 yu. Thus the 
angular distribution given in the table is for angles between 24.and 90°in the 
C-system, and the experimentally determined cross section is for the same 
range of angles. The minimum angle in the laboratory system is probably 
about the same in the cloud chamber experiment quoted above, and is 
much greater (30°) in the case of the counter experiment of Smith et al. 


in mb per steradian 


do 


da 
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We have made no correction for scanning efficiency because we believe 
that our scanning procedure has no appreciable loss factor. The results 
of Morris et al. are also not corrected for loss. 

The angular distribution which we have obtained is in good agreement, 
within the statistics, with that obtained by Morris et al. at an incident 
proton energy of 900 mev, and by Smith e¢ al. at 1000 Mev. It shows a 
strong forward peak in marked contrast to the isotropic distribution at 
400 mev (Hartzler and Siegel 1954). The deviations from isotropy begin 


Fig. 2 


437 MeV 


0 10 2Omeo0 405" 50. 60. 10 80 90 


Angle in C-system in degrees 
The angular distribution in the C-system for p-p elastic scattering at different 
energies. The full lines represent the results of Sutton e¢ al. at 437 Mev 
and Smith et al. at 590 and 1000 mev. The points marked A 
represent our results at 950 Mev ; the point -- is obtained from a 
combination of our results and those of Morris ef al. for the angular 
interval 0-20°. 
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to show at small angles at 428 mev (Hartzler and Siegel 1954) and at 437 
mev (Sutton et al, 1955). The work of Bogachev and Vzorov (1954) and 
of Selektor, Nikitin, Bogomolov and Zombkovsky (1954) in the U.S.S.R. 
(460-660 mey) and of Smith et al. (1955) in the U.S.A. (440-1000 Mev) 
shows that as the bombarding energy rises the differential cross section. 
increases markedly at small angles and decreases at large angles 
This is illustrated in fig. 2 in which we have plotted the data of Sutton et al. 
at 437 Mev and Smith et al. at 590 Mey and 1000 Mev, together with our 
experimental points. It is of particular interest to determine the behaviour 
of the differential cross section at the small angles. Because of the poor 
statistics in the range 0-20° (mean angle 14°) we have combined our 
results with those of Morris et al., normalized to our value at the cross 
section. We thus obtain the point shown which indicates that the 
differential cross section js still rising at small angles, although perhaps 
not as sharply as would be indicated by a direct extrapolation of the 
counter results of Smith et al. 

A study of the small angle (1°—5°) single scatterings observed in the 
course of the scanning indicates that the incident proton beam is not 
appreciably polarized. In any case our procedure has included the obser- 
vation of events in all azimuthal planes, and the results should therefore be 
independent of the polarization of the beam. 

Work is in progress on the analysis of the remaining 2p events in terms 
of (a) inelastic free p—p collisions and (b) elastic and inelastic collisions 
with bound nucleons. 
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ABSTRACT 

It has been found possible to measure the ionization of a cloud chamber 
track by photometry of the chamber photograph. The method is appli- 
cable to tracks, the photographic images of which consist of unresolved 
drop images. Its application does not interfere with accurate measure- 
ment of the magnetic rigidity of the particle causing the track. The use 
of the method has been investigated by studying tracks in a typical 
chamber and is illustrated by mass measurements on 3 K-mesons. The 
relative merits of this photometric method and the drop counting method 
of measuring ionization are discussed. 


§ 1. INTRODUCTION 

THE mass of a particle producing a track in a cloud chamber may be found 
by combining measurements of magnetic rigidity (which gives the 
momentum of the particle) and of ionization. The direct method of 
measuring ionization is to count the drops in a track. If this is to be done 
without sacrificing accuracy of momentum measurement, the tracks 
must be formed in a light gas and photographed at a comparatively large 
magnification (Cowan 1954, Blackett 1954). This has certain disadvan- 
tages which are discussed in § 7. 

The present experiment was therefore undertaken to see if it is possible 
to measure the ionization density of tracks using photographs in which the 
images of the individual drops in a track are not resolved. 

Three properties of the image of such a track might conceivably be used 
to measure ionization : the width, the structure and the mean blackness. 

Track width has been used (Burhop 1954) to measure ionization in high 
pressure chambers. It could not be used for the discontinuous tracks 
formed in a standard chamber. 

Using the structure of the track image, i.e. the variations in blackness 
along its length, one might hope to find a quantity depending on ionization 
but not on chamber and photographic conditions. Such a quantity could 
be used to measure the ionization density of a track without reference to 
that of other tracks. Preliminary studies in this direction were not 
promising and were therefore discontinued. 

The mean blackness of a track image depends on the track illumination, 
expansion conditions and film development as well as on ionization. In 
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general these change from expansion to expansion, and may be different 
for tracks appearing at the same expansion but in different parts of the 
chamber. Thus track image blackness cannot give an absolute measure of 
ionization. It has been found possible, however, to use it to compare the 
ionization densities of tracks formed and photographed under roughly the 
same conditions. 

The blackness of two track images may be compared in several ways, and 
a number of blackness-comparing parameters have been studied to see if 
they might be used to compare ionizations. It has been found most 
convenient to measure the blackness of a track image by a quantity which 
will be called the ‘ obscuration ’, 7’, of the image. This is defined by 


Hiei eani)it, nme ta Rs, He iC) 
where ¢ is the optical transmission of the track image and 1, is that of the 
portions of film on either side of the image. 


§ 2. APPARATUS AND EXPERIMENTAL DETAILS 


To test the feasibility of the method, tracks formed in the Pic-du-Midi 
chamber described by Armenteros e¢ al. (1951) have been studied. This is 
a cylindrical chamber, 28 cm in diameter, 10 cm deep, and is operated in a 
magnetic field of ca. 7500 gauss. The gas filling, which was at a pre- 
expansion pressure of 90 cm Hg, consisted of 85°% Argon and 15% Oxygen. 
Tracks were photographed on Ilford 5G.91 film at a magnification of 0-1. 

To measure the obscuration of a track image a microphotometer, using 
the Fiirth scanning principle (Firth 1943, Firth and Oliphant 1948), was 
constructed. A rectangular shaped light-spot was focused onto the film 
by a lens system attached to the prong of a tuning fork. The film was 
arranged so that the 60 c/s vibration of the fork caused the light-spot to 
scan a section of film normal to and including the track image. The scan 
amplitude was about ten times the width of a track image. Light 
transmitted by the film fell on a photomultiplier. Since this had a linear 
response, 

ie (Ged) LIV sin. Mthy oy oe El Paoaek m) 
where i is the multiplier current produced by light transmitted by the 
track image and i, is that produced by light transmitted by neighbouring 
portions of film. If R, and R are resistances carrying the multiplier 
currents, such that first i, and then i develop equal voltages, then 
T—R—R,/R. It was thus possible to make measurements of obscurations 
in terms of resistances only, i.e. independent of multiplier sensitivity, 
amplifier gain, incident light intensity ete. 

The length of the light-spot being 250, each measurement gave the 
obscuration of a 250 » long section of track image. A whole track image 
was measured in contiguous sections. A straight track image 2 cm long 
could be measured in 20min. The light-spot was 72 wide. (The 
track image width was such that 90%, of the drop images were contained 
within a width of 75 ww.) 
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§ 3, OBSCURATION AS A MEASURE OF IONIZATION 


If 7',/T',, the relative obscuration of two track images, is to be used 
for comparing ionizations, we must establish that this ratio is independent 
of the conditions under which the tracks are formed and photographed, 
provided that these are the same for the two tracks. 

For example, if the intensity of the light illuminating the two tracks 
is the same, 7',/7’, must not depend on this intensity. A fortunate 
peculiarity of the stereoscopic cameras used to photograph the Pic-du-Midi 
chamber allowed this point to be investigated. This is the fact that 
photographs taken by one of the cameras are slightly denser than those 
taken by the other. Thus each track is photographed twice at different 
exposures. This is similar to having the same track photographed twice 
under differing illuminations. 

In this way the effect of exposure on the image of a track of medium 
density has been compared with the effect on the image of a track of 
minimum ionization density. There is no evidence that the ratio 7',/T', 
for two such tracks depends on exposure. The measurements show, 
with 90°, certainty, that if a change of exposure alters 7', and 7’, by q%, 
any alteration of 7',/7', is less than 0-18 q%. 

There is, however, a more powerful test of the usefulness of 7',/7’, as 
a means of comparing ionizations. 


§ 4. Srupy OF PARTICLES OF KNOWN IONIZATION 


Suppose J, and J, are the ionization densities of two tracks formed 
and photographed under the same conditions. If 7',/7', is to be useful 
for comparing ionizations we must have 7,/7,=/(/, ; 1 2), where f(7,; I) 
is independent of chamber and photographic conditions. Of course 
obscuration can only directly measure the droplet density n of a track 
and not its ionization density. Now n=el, where « is the condensation 
efficiency, a quantity which will in general vary from expansion to 
expansion. Thus 7’,/7', can only be independent of chamber conditions 
if it is given by 7',/7,=f(1,/1.)=f(n,/n.). If this is to be true for all 
values of /, and J, the function must be a power law, Ue Deea(lailae 

The best way of deciding whether 7',/7', is independent of chamber 
conditions is the following. 7',/7', is measured for many pairs of tracks 
and a graph of 7',/7', vs I,/I, is drawn. Chamber and photographie 
conditions are the same for both tracks of a given pair, but they vary 
widely from pair to pair. If 7',/7’, depends only on J 1/1, then, despite 
the variation of conditions from pair to pair, all the points on the 
T/T’, vs I,/I, graph will lie on a single curve. In this case T/T. is a 
useful quantity for measuring relative ionization. If T',/T', depends on 
chamber conditions the points are spread over the graph. 

To provide data for such a graph (fig. 1), penetrating cosmic ray 
shower photographs have been selected, each containing an identifiable 
proton or deuteron track (identified by visual ionization estimation and 
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momentum measurement), together with at least one track at approxi- 
mately minimum ionization. lonizations were found by measuring the 
momenta and assuming the nature of the particles giving the tracks. 
_ Where possible a proton or deuteron has been compared with more than 
one minimum ionizing particle so as to reduce errors. 7, and J, are 
then weighted mean values for the minimum ionization tracks. 

No track contributing to the graph has an ionization greater than 
6 < minimum. 


Fig. 1 


The minimum ionizing particles have been carefully selected so that 
there is no ambiguity as to their nature. Protons have been excluded 
by using only negative particles. Supposed electrons have been used 
only if their momentum was less than 8 x 107 ev/c, at which momentum 
they are distinguishable from mesons by visual ionization estimation. 
Supposed mesons have been used only if their momentum exceeded 
5x10%ev/c. The number of electrons found in penetrating showers 
with momenta above 5 x 108 ev/c is very small (<1% of the total number 
of particles). All mesons have been assumed to be 7-mesons. For any 
that are really K-mesons the derived ionization is wrong. However, 
assuming that less than 20% of the mesons are K-mesons, it can be 
shown in every case that the chance that the derived ionization 1s wrong 
because the particle is really a K-meson is less than the chance that it 
is equally wrong because of measurement errors. 
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Only tracks on the same photograph have been compared. The spatial 
separation of tracks selected for comparison has been limited so that the 
effect on 7',/T7', of variations of conditions throughout the chamber is small 
compared with the effect of intrinsic blackness fluctuations (these are 
discussed below). 

When measuring the mean obscuration for a track, as when measuring 
ionization by drop-counting, it is statistically desirable to reject the 
contribution from the largest clusters. In drop-counting it is most 
important to decide on a limiting acceptable cluster size and to adhere 
to it ngidly. However, because of the saturation effects discussed below, 
the contribution of a large cluster to the mean obscuration is less 
important than its contribution to a drop count, and a rigid rejection 
criterion is not so important. Clusters have therefore been rejected 
simply on the strength of their visual appearance. The rejection criterion 
was chosen so that on the average one cluster was rejected for every 
5 cm of track. 

The points on this graph of 7',/7', vs 1,/I, are those shown in fig. 1 
with continuous error-bars. The errors in J,/7, derive from errors in 
measuring momenta. Errors shown vertically result from intrinsic 
ionization fluctuations which cause the mean blackness for any given 
track to deviate in general from the mean blackness for tracks of that 
velocity. The size of this error is found from the distribution about 
the mean obscuration of measurements for individual sections of a track. 
Compared with this error, errors arising in the microphotometer are 
negligible. All errors are given as standard deviations. 

The different points on the graph refer to tracks photographed under 
widely varying conditions. The image obscuration for a track of given 
ionization varies by an amount equivalent to an ionization ratio of 
>2-5 for tracks on the same photograph. Despite this wide variation 
in conditions, there is no evidence for a scattered distribution of points 
over the graph. They fit, within the errors, a single curve. Moreover, 
they are consistent with the power law relation between T/T, and 
1,/I, which is required if 7',/7', is to be independent of chamber conditions. 
Thus the facts that the points fit a power law and that no dependence 
of T',/7', on chamber conditions is apparent are in agreement. 

To test further if 7',/7', is independent of chamber and photographic 
conditions, it is useful to group together points which refer to tracks 
taken under similar conditions and see if there are systematic differences 
between the groups. Of course, with so few points on the graph the 
statistical weight of a conclusion based on grouping of points is never 
very high. Classifications have been made according to the image- 
blackness for a track of given ionization, according to the degree of film 
development and according to the degree of focusing of the track images. 
Jn no case is there evidence for systematic difference between groups. 


We have that 7,/7,=(I,/J,)". Making a least squares fit to the 
graph points n=0-73-+0-035. 
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In addition to protons and deuterons occurring in penetrating showers, 
7 positive secondary particles from A°-decays have been compared with 
approximately minimum ionizing particles. The latter were selected as 
_ described above. None of the A°-secondaries have ionizations greater 
than 6 x minimum. 

Results for the A°-secondaries can be treated in two alternative ways. 

Firstly, knowing the ionizations of the ‘minimum’ particles, we can 
deduce the ionizations of the A°-secondaries, using the 7',/7', vs I,/I., 
graph that has already been drawn. Combination with momentum 
measurements then allows the masses of the particles to be found. The 
measured masses, with standard deviation errors, are shown in fig. 2 (a). 


Fig. 2 


des 
a e 


They are all consistent with the mass of a proton. In figs. 2 (6) and 2 (c), 
the contributions to the total mass error from the errors in measuring 
ionization and momentum are shown separately. They are seen to be 
comparable and equal to ca, 200 m,. 

Alternatively, the A°-secondaries may be treated as known protons 
and a further set of points on the 7',/7', vs J,/I, graph obtained. If this 
is done the points shown by dotted lines in fig. 1 result. These new 
points are consistent with the old. A least squares fit to all the points 
gives n=0-72+0-028 in T',/T,=(L,/L.)”- 

We may summarize the above results as follows. 
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It appears that the relative image-obscuration of two tracks, 7',/T7's, 
may be used to measure the relative ionization density of the tracks. 
No evidence has been found for dependence of 7’,/7', on factors other than 
relative ionization density. 

We may therefore conclude that photometry is a valid and usable 
method for comparing the ionization density of tracks. Limitations to 
this conclusion are discussed below. 

It is very difficult to describe fully conditions obtaining in an unresolved 
track image and to calculate theoretically how its obscuration should 
vary with ionization. It has only been possible to calculate the behaviour 
of simple track models. Such calculations show that the observed 
obscuration—ionization relationship is reasonable. 


§ 5. RANGE AND ACCURACY OF PHOTOMETRIC 
JONIZATION MEASUREMENTS 


5.1. Range 


At ionizations of less than 6x minimum, 7',/7', for two tracks is 
independent of chamber and photographic conditions over a wide range 
of conditions. However, for the following reason, this cannot be true 
if 7’, refers to a track of high ionization density. 

There are a number of ‘ saturation ’ effects which tend to reduce the 
image-obscuration of a dense track: vapour depletion, photographic 
saturation, and the fact that obscuration, given by (1), approaches a 
limit of 1 for very dense tracks. If then a dense track is compared with 
one of low ionization density, 7/7’, will be reduced by these effects. 
The reduction will be greater if chamber and photographic conditions 
lead to dense tracks than if they lead to weak. Thus for two such 
tracks 7',/7', depends on chamber conditions. 

Even tracks of low ionization density contain clusters where saturation 
effects occur. The denser a track, the more dense clusters it contains. 
Thus 7',/7’, must depend to some degree on chamber conditions for all 
comparisons. The dependence for the present tracks is found to be 
unimportant if 7’, refers to an ionization of less than 6 minimum. 
It is found to become significant at ionizations 6x minimum. 

Thus, for the present photographs, photometric ionization measurement 
is essentially limited to ionizations of less than 6X minimum. 


5.2. Accuracy 
(4) Random Errors. As indicated in § 4, intrinsic fluctuations in 
ionization are the main source of error. Because of saturation effects 
for large clusters, the effect of fluctuations on a photometric measure- 
ment of ionization depends somewhat on chamber conditions. Under 
average conditions it is found that the error in a measurement of ionization 
is given by 


81/I=1-6/ni2 Pembina 6 fs setae (3) 
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where n is the number of primary ionizing collisions made by the particle 
in forming the considered length of track. Of course only relative ioni- 
zation is measured, and both the track of interest and those with which 
It is compared contribute to such an error. Thus, to measure the 
ionization of a track to the highest degree of accuracy, it will, in general, 
be necessary to use more than one comparison track. 

(6) Systematic Errors. In addition to random errors there is the 
possibility of systematic errors arising because conditions are not the 
same for compared tracks. 


Fig. 3 


Thus, fig. 3 shows (left-hand ordinate scale) that the obscuration for 
a track of given ionization depends markedly on its lateral location in 
the chamber. The right-hand scale shows the error introduced into a 
photometric comparison of the ionizations of laterally separated tracks. 

The points on fig. 3 were obtained by comparing separated tracks of 
equal ionization. Points shown with errors represent means of several 
such comparisons. Those shown without errors represent single compari- 
sons. 

The lateral effect is due primarily to changes in the angle through 
which the light reaching the camera is scattered at the drops. The 
intensity of scattered light depends markedly on scattering angle (Webb 
1935, Lander and Nielsen 1953) which differs by 20°. for tracks at the 
side and at the centre of the chamber. Effects of variations in light 
beam intensity are small compared with this effect. Figure 3 can be 
used to correct measurements of 7',/7', for laterally separatedt racks 
provided the corrections are not large. 3 

Spatial separations in other directions do not produce such marked 
effects. Thus, within the illuminated volume of the chamber, obscuration 
varies with the distance of the track behind the window of the chamber 
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by less than 10%. Nor is there a marked dependence of obscuration 
on the height of the track in the chamber. This shows that the vertical 
temperature gradient, maintained so as to give convective stability, does 
not markedly affect condensation efficiency. This has been confirmed 
by. Uttley (private communication) who has made drop counts in this 
chamber. The size of the temperature gradient is not known. For a 
larger chamber, a larger temperature difference between top and bottom 
would be needed to give the same convective stability, and vertical 
variations of obscuration, or droplet density in a drop counting experi- 
ment, would be more important. 

It is clearly convenient to be able to compare tracks on different 
photographs. 

Without standardized development, comparison of photographs on 
different lengths of film is impossible. For photographs on a single 
length of film it is a question of how steady chamber conditions can be 
kept. The charging time of the lamp condensers must be smaller than the 
least permitted interval between photographs, or the flash intensity 
depends on the actual interval. It is also essential to check the con- 
stancy of operating conditions, in particular the lamp voltage and chamber 
temperature, preferably recording these at each expansion. This was 
not done when the tracks studied in the present experiment were obtained, 
and so comparison of tracks in different photographs has in general not 
been attempted. However, measurements show that chamber conditions 
are frequently sufficiently constant to allow comparison of tracks formed 
at different expansions for periods of several hours. 


§ 6. PHOTOMETRIC MEASUREMENT OF MASSES 


To illustrate the use of the photometric method, results are given below 
of mass measurements on three K-mesons. The ionization of the particles 


Mass Measurement on K-mesons 


Photometric Measurement of Mass | Mass from visual 


Particle : : 
(in m,) with errors estimation 
| (a) (0) () (d) 
Charged V-particle 860 +175 +330 1580 + 160 
—200 —860 
Non-decaying par- 1000 +165 +290 850 + 250 
ticle from inter- —165 —320 - 
action 


Charged V-particle | 1220 at two standard deviations level 1660 + 160 
ee a eet aa ena re Sanne 


(a) Best value for mass. 
(L 
( 


») Error, putting 5p and 8/7 equal to one standard deviation (p=momentum), 
c) Error, putting 3p and 8/7 equal to two standard deviations. 

(d) dp put equal to one standard deviation: I has to be estimated subjec- 
tively, probably roughly equivalent to one standard deviation. 

The errors are combined quadratically. 
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was found photometrically in the same way as for the A°-secondaries. 
Momentum measurement then allowed masses to be determined. 

The table shows the results of the measurements. For comparison, 
_mass values based on visual ionization estimates are also given. 

All three particles are positively charged. For neither of the decaying 
particles is it possible to make measurements on the decay product. 

For a given chamber, the accuracy with which the mass of any particle 
can be found in this way may be calculated using (3), (4), and knowledge 
of the momentum errors in the chamber. 


§ 7. COMPARISON OF THE PHOTOMETRIC AND Drop Counting Mreruops 
OF MEASURING IONIZATION 


Both the photometric and drop counting methods of measuring 
ionization are much superior to the unreliable method of visual estimation. 

In both methods it is only possible to measure up to a limiting 
ionization ; in drop counting because drops become unresolved, in 
photometry because of saturation effects. Thus upper limit is raised if 
the photographic magnification is increased or the density of the chamber 
gas is reduced. The advantages of photometry over drop counting 
result from the fact that because resolved drop images are not required, 
either a smaller magnification or denser gas may be used. Comparison 
with drop counts made in the Pic-du-Midi chamber by Uttley (to be 
published) shows that for the same limiting ionization the number of 
drop images per cm? of track image can be 9 times as large in photometry 
as in drop counting. This means that either a denser gas giving 9 times 
as many ion pairs per cm or a magnification smaller by a factor 3 may be 
used. The advantages of doing one or other of these will be discussed 
in turn. 

7.1. Advantage of Using a Denser Gas 

In both photometry and drop counting 6//[=1-6/n" gives the error 
in an ionization measurement. Use of a denser gas reduces this error 
for a track of given length. The accuracy of momentum measurement 
on a track is probably not very dependent on gas density. 


7.2. Advantage of Using a Low Magnification 


Use of a low magnification makes photography more economical. Thus 
the area of film used in a single camera varies as m” (m—=magnification). 
Moreover, the camera aperture giving optimum depth of focus and 
resolution gives a depth of focus varying asm~*. Thus if several differently 
focused cameras are used to photograph a deep chamber, the total area 
of film used varies as m*. 

It is furthermore found, using the photometric method, that the 
relation between 7',/7', and J,/I, is little affected by slight defocusing 
of the track images. In the present experiment the method has been 
successfully applied to tracks varying in distance from the camera by 
2 to 3 times the depth of focus. 
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Other advantages of the photometric method are that it requires no 
modification of the chamber and can be applied to existing film records. 


7.3. Disadvantages of Photometry 


The first important disadvantage of the photometric method is its 
slowness. ‘This results from two causes: (1) the variation of conditions 
throughout the chamber must be investigated ; (2) not only the inter- 
esting track must be measured but also the comparison tracks, and this 
involves obscuration measurement, stereoscopic reprojection and momen- 
tum measurement. Thus it takes several hours to measure the ionization 
of a track photometrically. 

This disadvantage becomes less important as conditions are made 
more uniform throughout the chamber and steadier with time. 

A second disadvantage is that the requirements for the successful 
intercomparison of tracks formed at different expansions are somewhat 
more stringent in photometry than in drop counting. 

The comparison of the photometric and drop counting methods of 
measuring ionization may then be briefly summarized as follows. For the 
same degree of accuracy in ionization measurement, the photometric 
method is the less expensive. It is however the slower method. 
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ABSTRACT 


An alloy of composition close to Mn,AIC was prepared by induction 
melting. In the as-cast state the alloy was feebly magnetic at room 
temperature and the magnetic property was retained after homogenization 
at 1000°c. This alloy was composed of a single phase possessing an ordered 
face centred cubic structure with lattice parameter 3-869 A. The ordered 
structure is similar to that found in the alloy Mn,ZnC, the properties of 
which have been described elsewhere. The alloy Mn,AIC, is strongly 
magnetic at low temperatures having a saturation moment at 0°K of 
99-6 ergs/g/oersted equivalent to a Bohr magneton number 1-20 per Mn 
atom, the Curie temperature being 15°c. The paramagnetic susceptibility 
measured over the temperature range 65—-800°c, was found to vary with 
temperature in a manner different from a normal ferromagnetic material, 
the 1/x-7' graph possessing pronounced curvature concave to the 
temperature axis. 


§ 1. INTRODUCTION 


OVER a certain range of composition it is possible to obtain single phase 
ternary alloys of manganese, carbon, and aluminium possessing a face 
centred cubic structure showing pronounced ordering. The alloy of 
composition Mn,AIC exhibits the maximum ordering and has a detailed 
structure isomorphous with the alloy Mn,ZnC (Butters and Myers 1955). 
This latter alloy possesses interesting magnetic properties, being ferri- 
magnetic at low temperatures ; it is therefore of considerable interest to 
know the detailed magnetic properties of similar ternary alloys and this 
paper describes the properties of those alloys containing manganese 
aluminium and carbon. 


§ 2. PREPARATION OF THE ALLOYS 


The materials used in this work were manganese of 99-9°% purity 
donated by the Electromanganese Corporation of America, aluminium of 
99-99°% purity donated by the Aluminum Company of Canada and 
graphite of spectroscopic grade. 


*Communicated by W. Sucksmith, F.R.S, 
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The alloys were melted and chill cast under an argon atmosphere. 
Interest was centred on the composition MngpAl oC) but other alloys in 
which the carbon content was maintained at 20 atomic % and the 
proportions of manganese and aluminium varied, were prepared. In the 
as-cast state the micro-structure indicated pronounced coring ; this was 
removed by annealing at 1000°c for ten days, the homogenization taking 
place rather slowly. All specimens were sealed in evacuated quartz 
tubes for the homogenization treatment. 

The alloys are extremely unstable in air, particularly so under humid 
conditions, and it was essential that specimens be kept in a desiccator. 
The compositions referred to in later sections are the nominal compositions 
of the alloys; the actual compositions were determined by chemical 
analysis and found to be in good agreement with the nominal values, the 
results of the analyses are given in § 5 where the influence of the manganese 
in excess of 60 atomic °%, is discussed. 


§ 3. STRUCTURE AND LaTTicE PARAMETER 


Alloy specimens were prepared and examined under the microscope and 
by x-ray diffraction technique. The alloys were very hard and powders 
were easily prepared for Debye—Scherrer diffraction photographs. It was 
not necessary to anneal the powders to obtain sharp diffraction lines and 
all x-ray studies were made on unannealed powders. Fe Ka and Cr Ke 
radiations were used. Debye-—Scherrer photographs showed the alloy 
Mn,AIC to possess an ordered face centred cubic structure with lattice 
parameter 3-869 A. The intensities of the diffraction lines were measured 
using a Geiger counter spectrometer and were found consistent with an 
ordered structure similar to that already determined for the alloy Mn,ZnC 
(Butters and Myers, loc. cit.). In this arrangement the manganese atoms 
are located at the face centre positions, the aluminium at the cube 
corners and the carbon atoms at the body centre positions. 

The calculated intensities of diffraction lines arising from this structure 
and those actually observed for the alloy Mn,AIC are given in table 1. 
The temperature factor was omitted in these intensity calculations but the 
atomic scattering factors were corrected for any depression due to behaviour 
in the region of absorption edges (James 1948). 

The measured intensities given in table 1 may be compared with those 
found for the alloy Mn,ZnC. The superlattice lines are taken to be those 
with 2/?=1, 2, 5, 6, 9, 10, 13 and 14, and the geometrical factor occurring 
in the expression for the line intensity takes the values 


(f—fun—fc)? when Sh? odd 
and (f—funtfc)? when Sh? even 
where f is the atomic scattering factor for zinc or aluminium and fy, fo the 
same quantity for manganese and carbon respectively. 


Zinc possesses an atomic scattering factor equal, approximately, to the 
sum of the scattering factors of manganese and carbon, and in this system 


4 
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the intensities of the superlattice lines are controlled by the small contribu- 
tion of the carbon atom which leads to lines with Yh? even having greater 
intensity than those with Yh? odd. On the other hand the atomic scat- 
tering factor of aluminium is much less than the sum of the atomic 
scattering factors of manganese and carbon and superlattice lines with 
2/2 even have less intensity than those with Dh? odd. 

This reversal of the intensities of superlattice lines in the two alloys 
Mn;ZnC and Mn,AlC is to be expected on the basis of the assumed 
structure for these alloys and is actually observed in the Debye—Scherrer 
powder photographs but some of the diffraction lines arising from the 
aluminium alloy were too weak to be measured accurately with the 
Geiger counter spectrometer and a complete quantitative comparison is 
not possible. 


Table 1 


Line index Intensity 
Sh2 as he ern cern ek ng BN oe 
Measured Calculated 


not measurable 
26-5 

not measurable 

not measurable 
100 


1 
2 
3 
4 
5 
6 
8 
9 
0 
1 


— jl 


The measured and calculated intensities are based on a value 100 for the line 
ahe=11p 
Specimen composition Mngg.7Aly9.gC19.7 ; Lattice parameter 3-868 A. 
Chromium K« radiation was used for this intensity work. 


The occurrence of the intensity reversal between the lines with 21? 
odd and Xh? even is further evidence for the correctness of the assumed 
structure, being additional to the agreement between observed and 
calculated intensities found for the stronger diffraction lines of the two 
alloys. 

For a fixed carbon content of 20 atomic % the single face centred 
cubic structure occurs over the composition range Mn 60-69 atomic 
of Al 20-11 atomic%. The exact limits of this phase have not yet been 
determined, the primary interest being the magnetic properties of alloys 
with composition close to Mn,AIC. Preliminary investigations have 
shown however that for a fixed carbon content of 20% a two-phase alloy 
results when the manganese content is increased to 70%, or reduced to 
58% and it seems that although the Mn: Al atom ratio may be increased 


SE 2 


898 R. G. Butters and H. P. Myers on the 


to values considerably greater than 3, very little reduction, if any, can 
be made to this ratio and maintain a single-phase alloy. If, however, the 
Mn: Al atom ratio is kept constant at the value 3 a certain variation in 
carbon content is allowed. Single-phase alloys have been prepared with 
15%, carbon but carbon in excess of 20°% is dispersed as graphite. 
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The variation of the saturation magnetization with temperature for alloys of 
different composition. 


The alloys used for the magnetic work all contained a fixed carbon 
content of 20 atomic °%, and the composition was varied by increasing the 
manganese content. Very little change in lattice parameter was found 
for considerable: variation in composition which implies that in these 
alloys manganese and aluminium have almost the same atomic diameter 
of 2-73 A, the manganese being slightly greater in size than the aluminium 
atoms. 
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§ 4. Maagnetic PROPERTIES 


The saturation magnetization in an external field of 16 200 oersteds, 
and its variation with temperature were measured for alloy specimens 
using a Sucksmith ring balance. The paramagnetic behaviour above the 
Curie temperature was determined for the alloy Mn,AlIC using a similar 
but more sensitive ring balance. All magnetic measurements were made 
on coarsely powdered specimens which were found to saturate easily 
in the fields used. High temperature measurements were carried out 
im vacuo or a reduced atmosphere of argon. The use of coarse powders 


Fig. 2 
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The paramagnetic behaviour of Mn;AlC, y is referred to a gram of alloy. 
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required a special specimen carrier for the paramagnetic measurements ; 
this was made as an annular boat machined from graphite. Such a boat 
held about 120 mg material which was sufficient for the measurements. 
The carrier was diamagnetic and a correction was therefore determined for 
the temperature range used in these experiments (65-800°C). 
The variation of the saturation magnetization with temperature 1s 
shown in fig. 1, the values at 0°K being extrapolated value from Sms 
graphs and the Curie temperatures determined by plotting o? against 7’. 
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It is immediately seen that increasing the manganese content decreases 
the saturation magnetization and causes an increase in the Curie 
temperature. 

The variation of the paramagnetic susceptibility above the Curie 
temperature for the alloy Mn, AIC is shown in fig. 2, the shape of the graph 
is unlike that usually observed for ferromagnetic substances and possesses 
a curvature concave to the temperature axis similar to that found in the 
paramagnetic behaviour of the ferrite class of magnetic substances. The 
curvature of the 1/y, 7’ graph prevents a determination of the Curie 
constant for the alloy. 


§ 5. Discussion 


The alloy Mn,AIC possesses an ordered structure which is independent 
of heat treatment, the structure of the homogenized alloy being the same 
whether slowly cooled or quenched from 1000°c. It hardly seems correct 
therefore to class this observed ordering as typical of normal superlattice 
formation. Furthermore, the presence of partial ordering even in alloys 
differing widely from the optimum composition Mn,AIC suggests that 
another mechanism must be found for the observed regular structures. 

For alloys containing a fixed carbon content of 20 atomic °%, the 
manganese aluminium atom ratio may be increased to values considerably 
greater than 3, the optimum value, but a second phase is produced when 
this ratio is decreased. This implies that the aluminium atoms cannot 
occupy face centre positions in this structure. A possible interpretation 
is that the carbon atom prefers manganese nearest neighbours and a 
consideration of atomic diameters leads to the conclusion that the carbon 
must be considered bonded to the manganese atoms. (The behaviour 
of the Mn, Zn, C alloys with the same structure is similar in this respect.) 
The regularity of the structure may be understood if it is assumed that 
vacant d orbitals of the manganese are used in the bonds with the carbon 
atoms ; increasing the manganese content would not affect this mechanism - 
but the low manganese limitof the phase can be appreciated, for if alumin- 
ium occupies a face centre position there is no simple way whereby the 
aluminium can provide d orbitals to form bonds with carbon. 

If then the tendency to form an ordered structure is governed by bond 
formation the preservation of the order after different heat treatments 
can be understood. On the other hand, the equivalence of the six 
manganese atoms in the unit cell of Mn,AIC demands that the carbon 
on the average be bonded equally to each. We do not wish to discuss a 
detailed model for this process since we believe that additional data 
obtained from neutron diffraction experiments will enable a more precise 
appreciation to be made. 

The saturation magnetization varies with temperature in a normal 
manner and there are no indications of any anomalous behaviour as found 
in Mn;ZnC. The Bohr magneton value of 1-21 is similar in magnitude 
to the value found in the zinc alloy and very different from the value 4 
found for manganese ir, the Heusler alloys. The paramagnetism of 
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Mn,AIC above the Curie temperature is definitely unusual for a metal 
and indiates ferrimagnetic rather than ferromagnetic properties. For an 
arbitrary value of Curie constant the equation derived by Néel (1948) 
_for the paramagnetic susceptibility of a ferrite may be fitted accurately 
to the experimental points. If, however, ferrimagnetism is a proper 
interpretation of the paramagnetic behaviour then there is an unusually 
large difference between the ferromagnetic Curie temperature and the 
6, of Néel’s treatment (0, will be ca. 160°c) this difference being about 
150° compared with the usual 40 to 50° for this difference. 

Other than the low Bohr magneton value for the manganese atoms and 
the form of the 1/,—7' curve there is no direct evidence for a ferrimagnetic 
interpretation but the manganese-manganese distances in the alloy 
Mn,AIC should favour negative interactions. In Mn,AlC and Mn,ZnC 
these distances are 2-73 A and 2-76 A respectively, and at low temperatures 
Mn,ZnC is known to be ferrimagnetic. Furthermore, Guillaud (1953) 
using data obtained from Mn,N, MnSb, MnAs, and other compounds 
suggests that negative interactions will be found when manganese atoms 
are separated by less than 2-814. An interpretation of the results in 
terms of ferrimagnetism would therefore seem appropriate. 


Table 2 


nae en eee 


Be cition MnggAlogCoo MneteAle Oss) = | MasAls0= MnypAl,gCop* 


feasured alloy 
ne Ii oI AO. Min see Aleeee Conn Miele ay Cone — 
composition 60°7 19°6~19°7 61°95 17°97 20°08 64°05 15°61 ~ 20°34 


— 


F This alloy contained a small amount of a second phase, the manganese contributing to the 
netic moment would therefore be less than 70 atomic °/ but can not be estimated. accurately 
Increasing the manganese content beyond 60 atomic % causes an 
increase in the Curie temperature and a decrease in the saturation moment 
of the alloy. Assuming that only the manganese atoms contribute to 
the magnetization the Bohr magneton number (pz) per manganese atom 
was calculated for each alloy. For this purpose the alloys were analysed 
chemically for each constituent element. The totals in all cases were 
within the limits 100-0-5% by weight and the atomic percentages were 
determined by scaling to 100%. The nominal and determined compo- 
sitions with the calculated p,; values are given in table 2. The variation 
of pz, number with the manganese content in excess of 60 atomic % is 
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shown in fig. 3. Since the alloy containing 70 atomic % Mn had a small 
amount of a second phase present the values are not sufficiently accurate 
to be compared with the other alloys. The remaining alloys give results 
which show a linear decrease in moment with manganese in excess of 


Fig. 3 
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The variation in Bohr magneton value with manganese content in excess of 
60 atomic %. 


60 atomic °% and the gradient of the graph indicates that the magnetiza- 
tion is reduced by 5--0-5 pg per manganese atom in excess of 60 atomic %. 
Manganese atoms in excess of 60 atomic °%, must replace Al atoms at cube 
corners and the magnetic results indicate that they probably exist in 
3d,4s, states and give rise to a moment antiparallel to that of the 


manganese atoms at face centre positions. 
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SUMMARY 


Pople’s recent criticism of the results obtained, for the high temperature 
susceptibility of a permanent dipole lattice, by Van Vleck and by 
Rosenberg and Lax is shown to be unjustified. Two alternative deriva- 
tions are given of the relations obtained by Van Vleck e¢ al. It is shown 
that Kirkwood’s correlation parameter, g, depends upon Pures 
interactions. 


§ 1. INTRODUCTION 


THE problem of the electric polarization of a system of permanent dipoles 
arranged on a lattice has received considerable attention. Van Vleck 
(1937, 1940) has given an exact classical treatment of the problem and 
has obtained an expansion for the high temperature susceptibility in 
powers of (1/7’) as far as the terms in (1/7')?. Rosenberg and Lax (1953) 
have reorganized Van Vleck’s treatment and have extended his expansion 
to include the term in (1/7’)*. Recently Pople (1953) has stated that the 
results obtained by Van Vleck and by Rosenberg and Lax are in error. 
In the treatment given here it is shown that Pople’s criticism is unjustified, 
and that the results obtained by the previously mentioned authors 
are correct. 
§ 2. FORMULATION OF THE PROBLEM 


The starting point of this discussion is the general relation for the 
dielectric constant, «, of a uniformly polarized body, of volume v, 


(e—1)=(4nP/E)=(4n(M.e)oH)t . . . . . (2.1) 


where P is the dipole moment per unit volume, E (~e/) is the uniform 
potential gradient existing in the volume v, and M is the total instantaneous 
dipole moment of the volume. 

For simplicity it will be assumed in the following that the volume, », 
is spherical in shape. This spherical specimen of dielectric may be in 
vacuo or surrounded by an infinite volume of its own medium. It will 
be shown that identical expressions for (e—1) are obtained in both cases. 


* Communicated by the Author. 
+¢ ) denotes an average. { )° indicates that the average must be evaluated 


when the external field is zero. 
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Provided the spherical volume is sufficiently large to be considered as a 
macroscopic specimen, we may obtain the following relations from (2.1), 
thus : 
(a) sphere in vacuo, 
3[(e—1)/(e+2)]=[47(M .e),a,/0F] . . . - (2.2) 
where & denotes the external field in which the sphere is situated, and 
(b) sphere embedded in its own medium, 
(e—1)=[47(M Se /ebe .) e e 
Frohlich (1949), following the work of Kirkwood (1939), has shown that 
(2.2) and (2.3) may be written, ignoring saturation, in the following form, 
(a) B[(e—1)/(e+2)]=[4r(M2)o,/30kT] . . . (2.4) 
(b) [(e—1)(2e+1)/3e]=[47(M2)°/30kT], . . . . (2.5) 
where k is Boltzmann’s constant and 7’ the absolute temperature. 
(M2)°,, is the mean square fluctuation of the dipole moment of a 
dielectric sphere in vacuo and <M?)° is the mean square fluctuation of 
the dipole moment of a dielectric sphere embedded in an infinite volume 
of its own medium*. It is particularly to be noted that (2.5) has been 


obtained by considering the surroundings of the sphere on a strictly 
macroscopic basis. 


The problem is therefore, to evaluate (M?);,,. and (M?)°, and to show 
that the values so obtained, when inserted in (2.4) and (2.5) yield the same 
result for (e—1). Attention will be confined to lattices of tetrahedral 
symmetry. 

§3. EVALUATION OF (M?)°,. 


The evaluation of (M*)%,., will be first discussed as it is somewhat more 
straightforward than that of (M*)°. The procedure and notation adopted 
here are the same as in Pople’s (1953) work. We may write 


N 
(M? vac (2 Be Ue; aoe Ne ONS ies ext (3.1) 
U9 


where N is the number of dipoles, all of equal magnitude, and p is the 
dipole moment of a molecule. 


We assume that all dipoles in the sphere are equivalent. Consequently, 
N 
(M? yvac= NV C2". ui + Hi Dyno" soe ee (3.2) 
] 


To obtain an expression for (M?)\,, we must evaluate the following 
integral 


-N , 
(M2 yiuo=| [2 - UB; exp (—Ujk) 42 | | exp (— U/kT’) 1a | é (3.3) 


where U is the total interaction energy of the system, for a given configura- 
tion, and dQ is an element of configuration space. In the case which 


i ee 


* Mean square fluctuation since (M )vac=<{M )°=0. 
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we are considering, U takes the simple form 


N 
DEA pee meee we. 3s (3.4) 
; ij 
- where, in Pople’s notation, 
W p= op Mins Hsp a Mee mes. fs. «(o.0) 


<\)) is the dipole interaction tensor, 


Cee ton, —3R, E41. oe lo nome (3:0) 


tIh 


R;; is the radius vector joining lattice site i to lattice site j7. Greek 
suffixes denote Cartesian components of a tensor and may take the values 
1, 2 or 3; repeated suffixes are to be summed. 

After expansion of the exponentials in the integrals in (3.3) and using 
the facts that 


| ont, 005 tor 147, 


and 
| -tee-ty.wydQ=0 for all i, j,k, and], 
we obtain 
NN 
(M2 Nac =NueLN | 53 See 10 E j a 
HAY) 
+terms involving (1/7')? and higher powers. . — se (sal) 


Note that only finite sums are involved in (3.7), in contrast to Pople’s 
calculation, and therefore the order of summation is immaterial. 
Integration of (3.7) over 2 yields 


N N 
(M? onc Ne? + (1/27) (Np /keT?] ras 2 eo ecb an a( 3.8) 
Mea 
where 
ee oh patltaae(s(Rip Ray) ty |, eo) 


Equation (3.8) may be written 
(M22, = Np2-+(1/27)[ Nu (R22?) | 22 LEP _Tepep |. . 8.10 


But because the lattice is assumed to have tetrahedral symmetry and 
because the specimen is spherical in shape, 


N 
Ee ON eee (3.11) 
ij 
and so (3.10) becomes 
(M? vac Nu? — (2/9) es = Fa joe eee (3.12) 


If the substitution b=[47Np?/3vkT'] is made, (2. * C0 in conjunction 


with (3.12), 
$[(e—1/(e+-2)|=p— By? ee nea, (a.k3) 
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N 
Po BE Vvj4aNY SR one 
k( #4) 

is a number depending on the type of lattice. For the simple, body- 
centred, and face-centred cubic lattices B has the values 0-10642, 0-09196 
and 0-09154, respectively, for large NV (Pople 1953). Equation (3.13) is in 
agreement with the result obtained by Van Vleck (1937)* and Rosenberg 
and Lax (1953). Equation (3.13) may also be written in the form 


(e—l)= bf 145 ae Bye aan al &{3.105) 


§4. EVALUATION OF (M?)° 


The evaluation of (M?®)° differs from that of (M?)°,, because one has 
to include in the potential energy term of the Boltzmann integral the 
interaction energy of the sphere with its surroundings. This energy of 
interaction, W, will be computed by considering the surroundings of the 
sphere on a macroscopic basis. It will be assumed, as Frohlich (1949) 
has done, that W and M are related in the following manner, 

W =— (4)M?[477/3v][2(e—1)/(2e+1)]=—pM?. el) 
Equation (4.1) may be derived on the assumption that M may be treated 
as a point dipole situated at the centre of the spherical region. This 
last assumption does not give rise to any error in the terms of the 
expansion for (e—1) evaluated here. Further, the use of averaged 
quantities in the expression for W may be justified. Proofs of these 
statements will, it is hoped, be given elsewhere. 

Thus, for the problem in hand, U is given by 


N 
= Wee aige oc. 2 gk SO ae a ae 

i>j 
Pople (1953) attempted to evaluate <«(M?)°, but used an incorrect 


expression for U. The value of U given by (4.2) contains, besides the 
N 

sum 2'w,,, which in this calculation is finite, the factor W, which was 
i>j 

omitted in Pople’s work. 


From (4.2) we have 


 N N 
(M?)°—= NV | | Ey + by exp = W+ 5 w,,) nev 10 
j 


i>j 
re N aa 
x | Jexp{—ar+ 2: Wj) ur} ao | a eet) 
i>j 
A slight digression must be made here to consider the temperature 
dependence of p (see eqn. (4.1)), because we are only interested in those 
terms of the expansion of (M®)° which do not involve powers of (1/7) 


* The expansion given by Van Vleck (1940) for (e—1) appears to contain a 


misprint. Equation (3.15) may be derived directly from the result given in his 
earlier paper (1937). 
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greater than the second. Clearly, the leading term in (4.3) will be Np, 
and so (2.5) must be of the form 


[(e—1(2e+1)/SeJ=Y[1+A/T)+o(/TP+...]; 9. . (44) 
by suitable manipulation of (4.4) it may be shown that 
p=%[1-+ terms involving powers of (1/T)]. . . (4.5) 
A calculation, similar to that given in § 3, yields the result 
(M2)°=Np2-+ (2/9)[NA4/kT I[4n/30—By®. . (4.6) 
And so we obtain 
[(e—1)(2e+-1)/3e]= [y+ (2/99 BY +...J. 0. (4.7) 
Equation (4.7) may also be written, 
2 
(—1)=y| 1 a ia B+... |): eee (458) 


This result is the same as (3.15) which, as previously stated, is in 
agreement with Van Vleck (1937) and Rosenberg and Lax (1953). The 
expression derived, by Pople for <(M?)° is in fact equal to (M?)°... 
Replacing «M?)° by (M?)°,, in (2.5) gives Pople’s expression for (e—1), 


§5. Discussion 


In correspondence Pople has raised certain objections to the method 
of calculation employed in this paper. He states that the relation 


(th; + Uy ao BE Ye; B; Exp (— U/L) a2 | | exp (—U/kT) 10 ah 
(3.3a) 


is strictly valid only when NV, the number of dipoles, is very large and 
that (3.3) will be valid only if both lattice site 7 and lattice site 7 are 
surrounded by a very large number of dipoles. He considers that one 
should write 


dir Deao= | J-te By exp (—U/kE) 12 || | exp (—ujeryae |" +x 
(5.1) 


where X > 0 as N + o, and that it must be proved that NX is negligible 
in the present calculation, where V, although large, is finite. 

It would seem that the following discussion overcomes these objections. 
For a system of n particles whose total energy, H, is constant, and where 
the number of particles with energy «,; is n,;, the probability P of finding 
the system in a particular state is given by 


P=(n!)(n71 2,! nz! ...2,!.-.)7 1. const. ee 0:2) 
for Maxwell—Boltzmann statistics. In the usual way (Tolman 1938) one 
finds that the most probable value of 1, is given by 

n,=exp (—a—fe,), Meee ues (Ord) 
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where « and f are constants, on the assumption that ; 
log n,; !=n, log n,—nj. eM ie Be tad ied, 
Equation (5.4) will only be a valid approximation when ”, is very large. 
However, one may use a better approximation than (5.4) for log n,!, 
namely 
log n,!=n, log n;—n,+(4) log 27n,. ee ee (E59; 
To this approximation one finds that the most probable value of x, 
becomes 
n,=exp (—a—Be,) exp [—(—«)—(n—f)e,—E/n,], ~ (5.6) 
where ¢, 7 and € are constants. In the limit of large 7,, 
C>a,»—>B8 and é/n, > 0. 
Equations (5.3) and (5.4) may be rewritten in the following manner : 
n,=exp (—a—fe,)=fo 53 0 gt peghne Ma Sen ee 
nj=fy exp [—(S—«)—(n—B)e;— /n J=f. - . (5.8) 
The error term, X, in (5.1) referred to by Pople, corresponds to 


fofexp [—(¢—«)—(n—B)e,—E/n,]—]}. 
Pople’s objection may be stated thus: in the calculation of <p; -; Vac 
(§ 3), fo has been used to determine the distribution function, and the 
factor (eqn. (5.8)) 
exp [—(¢—a)—(y—B)e,—€/n,] 

has been assumed to be equal to unity. However, if the error in setting 
this factor equal to unity were appreciable we would then expect that 
the quantity obtained for (uw; .w;)f,o would depend upon the value of N, 
the total number of particles of the system. It will be seen that, in the 
expansion obtained for (u;.u;)fac, the only term which is affected by 


N 

the magnitude of N is that involving the summation » R,,~®. It is 
k( At) 

of course necessary that N should be sufficiently large so that this sum 


N 

should approach its limiting value. 2 R,,~° will not have exactly the 
k( 4%) 

same value at every lattice point, but, because of the rapid convergence, 

the value of the summation will be very nearly equal for all lattice sites 


except those near the surface of the specimen. The error arising from the 
N 


assumption that 2 R,,~*° is the same for all lattice sites of the specimen 
ea 
may be made Panialiael small by taking NV sufficiently large. 

By similar reasoning it would seem that Pople’s second objection can 
also be answered. It will be seen in §3 that (u,.u;)°,., only depends 
upon £;;, that is on the relative positions of pw; and p,, and not on their 
positions relative to the surface of the specimen. This result arises from 


the symmetry of the lattice and the symmetry of the shape of the 
specimen, 
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The fact that the expansion for (e—1) does not depend upon NV implies 
that the dielectric constant of a macroscopic body does not depend on its 
volume. Further, the fact that (u,.u,)°,. depends only on the relative 
positions of p; and w, and is very nearly the same for all i and j means 
_ that the body may be considered as homogeneous from a macroscopic 
point of view. In their paper, Rosenberg and Lax (1953) have shown 
that the factors in the expansion for (e—1) which depend upon shape 
cancel out explicitly. 

Therefore it would seem that (3.3) is valid for macroscopic bodies. 


$6. CONCLUSIONS 


It is concluded that Pople’s attempt to calculate (M2?)° was unsuccessful, 
and that the result he obtained was in fact equal to (M2)°,,. He ignored 
the fact that the relation (2.5) (eqn. (1.2) in Pople’s paper) was obtained 
on the basis that the surroundings of the sphere are to be treated as a 
continuum of dielectric constant «. 


The relation 
[(e—1)(2e+1)/3eJ= PLL + (2/9)2—Bye+ J. (4.7) 


to the approximation considered in this paper, differs from Onsager’s 
relation, 


[(e—1)(2e+1)/8eJ=h Cate ee on Os) 
(Onsager 1936) by the terms [(2/9)4?— Bys?]. 
Kirkwood (1939) has developed a theory of the dielectric constant 
which yields the formula 


[(e—1)(2e+-1)/8eJ=gs Ae Dade Wie Sern (Or) 


where g, known as the correlation parameter, is a factor which takes into 
account the effect of short-range non-dipolar forces. Kirkwood expressed 
the opinion that g would be equal to unity for pure dipole-dipole 
interaction. But, although only dipole-dipole interaction has been 
considered in this paper, one finds that 


FEAT IE ay eee (2) 


Thus g may differ from unity even in the absence of non-dipolar short- 
range interaction. Attention has already been drawn to this point 
(Frenkel 1946), nevertheless it seems desirable to emphasize the fact 
that g does depend upon dipole-dipole interaction. 

Finally, it must be stated that the conditions necessary for the 
convergence of (3.15) are not known, but presumably the expansion will 
hold for sufficiently small values of ¢. 
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ABSTRACT 


Systematic measurements of thermoelectric force and _ electrical 
resistivity of dilute alloys of iron in copper have been made at close 
intervals of temperature between 4:2°k and ~ 50°K. These show that 
iron gives rise to particularly large anomalous thermoelectric effects 
and minima in the resistance when dissolved in copper. The anomalous 
resistive component calculated by methods described in this paper does 
not follow an exponential or power law dependence on temperature, but, 
within certain limits, follows an empirical law dependent on a character- 
istic temperature. Taken at its face value this would suggest that the 
resistance minimum, in this case at least, is the result of some transition 
phenomenon. 

The size of the resistance minimum in these alloys is notable; for 
instance, the anomalous increase in resistivity between 27° and 4-2°K in 
an alloy containing 0-044 atomic °% Fe is equal to the thermal scattering 
component found in pure copper at about 50°K. 


§ 1. INTRODUCTION 


Iv has been known for some time (Borelius, Keesom, Johansson, Linde 
1930, 1932) and confirmed recently by us that as little as 0-004 atomic % 
Fe in solid solution profoundly modifies the thermoelectric power of 
copper at low temperatures. We have also found that this quantity of 
iron modifies the electrical resistivity of copper at low temperatures, 
giving rise to aresistance minimum of considerable size. It was, therefore, 
of interest to continue our work on the electrical properties of dilute 
copper alloys (MacDonald 1952, MacDonald and Pearson 1952, 1953, 
1954 a, 1955 a, b—I and II of the present series, Pearson 1954, 1955) 
by examining in more detail the behaviour of very dilute alloys of iron 
in copper. 

The equilibrium solid solubility of iron in copper decreases very rapidly 
with temperature, being ~ 0-6%, by weight at 800°C, 0:1% at 600°C and 
negligibly small at room temperature (Bauer and Hansen 1934). In 
order to avoid the possibility of effects due to coherently precipitated 
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iron (Reekie, Hutchison and Hetherington 1953) we have examined only 
alloys which appear to lie in the homogeneous solid solution range at 
the annealing temperature. It is, however, interesting to note that 
iron, lead and bismuth, all of which have a rather small solid solubility 
in copper, show effects on the electrical properties of copper qualitatively 
somewhat different to those of Ga, In, Si, Ge and Sn, where the solid 
solubility is much larger (greater than 5 atomic % at the annealing 
temperature of 530°C). 


§2. ExPERIMENTAL METHOD 


As in earlier work we have prepared alloys from a source of copper 
known as ‘JTH’, which had the lowest residual resistance of any 
polycrystalline copper which we had examined. This copper contained 
~0:002 % by weight Fe (present as an oxide), some oxygen which 
exceeded the solubility limit (~ 0-0015% at 500°c), and traces of Ca. 

A range of alloys was prepared from a basic alloy containing ~ 6% Fe 
and JTH copper by melting in alumina crucibles under hydrogen which 
was pumped off before the alloys solidified. In this way, the oxygen 
present was removed and all of the added iron entered into solid solution 
in the copper. The alloys were cold-worked and homogenized at 1050°c 
to 1065°c for a few hours or more. After cold-rolling and cutting into 
strips 2mm x 0-1 mm, the specimens were annealed for a minimum of 
14 hours at 530°c and quenched. The electrical properties were measured 
by methods previously described (MacDonald 1947, MacDonald and 
Pearson 1953, 1954 b). 

The iron content of the alloys has been determined by chemical 
analysis and checked by spectroscopic analysis but, as we have already 
pointed out elsewhere, when very small additions of specific solutes are 
made the most satisfactory method of determining the amount in solid 
solution is by measurement of the residual resistivity of the alloy. 


§3. MEASUREMENTS OF THERMOELECTRIC FORCE AND 
ELECTRICAL RESISTIVITY 


The absolute thermoelectric power (S=dH/dT’) of pure copper is very 
small up to 10°K after which it grows with positive sign. We find that 
even the smallest additions of iron (~ 0-002 wt. %) make the absolute 
thermoelectric power negative, increasing in magnitude very rapidly 
with temperature at first, becoming constant in the region of 20°K and 
at high temperatures decreasing again slowly in magnitude. This 
agrees with the findings of Borelius e¢ al. (1932), although at 15°x, say, 
where S is relatively constant for a given alloy and thermal scattering 
is still small, our values of S,,., are a little larger. Figure 1 shows how 
Si59x Increases to a maximum at 0-05 to 0-07 atomic % Fe and then 
decreases again slowly. As the residual resistance of alloys increases 
linearly with added iron content up to ~ 0-05 at. % Fe (fig. 2), the 
flattening-off of the S,;., versus 7,,, curve before this concentration 
cannot be the result of precipitation from solid solution in our alloys 
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annealed at 530°c, and a fortiori it cannot be the case with the alloys 
investigated by Borelius et al. which were quenched from 800°c, 
Fig. 1 


K INALV/°K 


Sis° 


0 0.1 0.2 0.3 0.4 0.5 
) AT. % IRON 
Variation of absolute thermoelectric power at 15°K with iron content of 
dilute Cu—Fe alloys. 
V present work ; alloys quenched from 530°C. 
© Borelius et al. (1932) ; alloys quenched from 800°c. 


Figure 2 shows the variation of residual resistance with analysed iron 
content of the alloys. The analyses of these very dilute alloys are 


io, 2 
Fig. 2 
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Variation of residual resistance, /yes, with analysed iron content of alloys. The 
analyses are reliable generally to +20%. 


uncertain to 120%. We find that the increase of residual resistivity of 
copper per atomic %, of Fe dissolved in solid solution is ~ 10-4 x 10~° ohm 
em which compares with the value of 9-3x10~§ ohm cm obtained by 


392 
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Linde (1932) from measurements of the resistivity of alloys down to 
78°x. Iron, therefore, has a scattering cross section of 1-8 x 10-79 em? 
in copper, which, it should be noticed, is a rather large scattermg cross 


section. 


Fig. 3 
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Comparison of the behaviour of dilute Cu-Sn and Cu—Fe alloys. 

Top: | Variation with residual resistance of absolute thermoelectric 
power at 15°K,. 

Middle ; Variation with residual resistance of relative rate of resistance 
increase below the resistance minimum. 

Bottom : Variation with residual resistance of the ‘ absolute’ rate of 
resistance increase below the resistance minimum. 

We have taken the slope, dR/d7’, somewhat above 4°k, where it is 

relatively constant, as a measure of the rate of increase. This parameter 

may also be regarded as a measure of the ‘size’ of the resistance 

minimum. 
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The size of the resistance minimum found in these alloys is particularly 
notable. White (1955) has examined the resistivity of the alloy shown 
in fig. 5 (Cu+-0-044 atomic % Fe) between 4:2 and ~ 0-1°K. Combining 
these measurements with ours, we find that the anomalous increase in 
resistivity between 7’',,,, and 0-1°K is greater than 16% of the resistivity 
of the alloy at Tin. 

In order to compare the behaviour of iron as solute in copper with 
that, say, of tin, we show in fig. 3 the absolute thermoelectric power at 
15°K, the relative value of the resistance minimum and the ‘ absolute ’ 
size of the resistance minimum as functions of the residual resistance 
(solute concentration) of the alloy.* The absence of the sharp cusp 
such as is found at 0-005 to 0-01 wt. % solute for alloys of Sn, Ge, In 
and Ga in copper is notable. The data shown in fig. 3 were obtained from 
measurements of thermoelectric force and resistivity of the alloys 
continuously between 4-2 and ~ 50°K. 


$4. Discussion 


In considering the resistance minimum in these alloys, our best approach 
seems to be to separate the total resistivity into four parts, such that 
Tootal="r+rpt+rotrs, where rp is the thermal scattering which would 
presumably be the only resistance in an ideally pure copper sample, 
rp is scattering due to various physical and residual chemical imperfec- 
tions, 7g is the normal scattering due to the specific added chemical 
component causing the minimum and r, is the anomalous scattering due 
to this component. rp and rg are usually considered together as 7inpurity 
and are presumed to be temperature-independent. They may, however, 
be separated if 7p is taken as the residual resistivity of pure copper which 
does not show a resistance minimum above 4°k. We may obtain 
(rp+1c)+r4 for a copper—iron alloy by subtracting rp as determined 
from pure copper. 

In figs. 4 and 5, which show the resistance-temperature variation of 
two Cu-Fe alloys, rp at any temperature is given by ordinates within 
the region ABC and values of (rp+rg)+7,4 are given. by the curve DEF.7 


* We have actually taken the resistance gradient a little above 4°K (in 
which region it is relatively constant for a particular alloy) as a measure of the 
size of the anomalous rise of resistance, as this parameter does not appear to be 
directly influenced by variation of the position of the minimum (7Z"njin). 

+ We assume that Matthiessen’s rule would hold for copper-iron alloys in 
relation to pure copper if there were no anomalous scattering component 
present at low temperatures. Slight variations do in fact occur which are all 
in the same sense relative to pure Cu. Similar deviations have previously been 
reported for copper alloys by van den Berg (1938) and by Griineisen (1933). 
These variations have the effect in general of making the curve for the anomalous 
component DE (figs. 4 and 5) apparently rise again beyond E. They will also 
affect the apparent position of E but we shall not consider this question further 
here as any changes produced generally lie within the overall accuracies attri- 
buted to the measurements. 
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As (rp-+rg) is postulated as temperature independent it is reasonable to 
assume that the anomalous component 7, is given at any temperature 
by ordinates within the region DEG, ie. r 4=" total —'2— (Te +1¢)- 
Examining the temperature dependence of the anomalous component 14 
we find that it does not follow a simple exponential or power law of the 
type r4=f1(rc) exp (W4/RT) or r4=fo(r¢)(0/T)” (fig. 6). It does, however, 
appear to be represented surprisingly well by the relationship 
ra=f(ro)(O—-T)” ; 
between 6 and 4-2°x* (fig. 6) where f(r;) is constant in any particular 


Fig. 4 
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The resistance minimum in an alloy containing ~ 0-044 atomic %/ Fe. 
©, Curve DH: The resistance of this alloy as a function of temperature. 
v, Curve BA: The resistance of two‘ pure ’ copper samples as a function 
of temperature. 


A, Curve DE: Variation with temperature of the anomalous resistive 
component, as deduced in the text. 


* It appears that this relationship is not obeyed at lower temperatures since 
White’s (/oc. cit.) results show a marked flattening of the resistance below ~2°x, 
becoming practically temperature-independent below 1°x. 
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alloy, 7’ is the absolute temperature, @ is the temperature at which r 
apparently becomes zero (point H, figs. 4 and 5), @ being a function nf 
the iron content of the alloys, and n is approximately 2 for alloys 
containing between ~ 0-005 and ~ 0-05 atomic % Fe. Although it 
would be of interest to proceed further by examining in some detail the 
variation of @ and f(rg) with r¢, we do not believe the accuracy of the 
present data justifies this step. f(r¢), however, does appear to be more 
or less a linear function of rg. 
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The resistance minimum in an alloy containing ~ 0-005 atomic % Fe. 
©, CurveDH: The resistance of this alloy as a function of temperature. 
vy, CurveBA: The resistance of two ‘ pure ’ copper samples as a function 
of temperature. 
A, CurveDE: Variation with temperature of the anomalous resistive 
component, as deduced in the text. 


We can provide no specific theoretical interpretation for such a relation- 
ship at present. However, if one were to accept the expression at its 
face value, then we should be forced to recognize the anomalous resistance 
as some kind of transition phenomenon setting in at a definite tempera- 
ture, @. Established phenomenon of this kind (e.g. ferromagnetism, 
superconductivity, the liquid helium I-I transition) always involve 
some cooperative behaviour in the assembly, and the * phase’ changes 
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are generally regarded as being of the second order. The thermoelectric 
power, because of the Thomson relationships, should be a particularly 
valuable parameter for examining further the validity of such an inter- 
pretation. For a second order process one might expect a discontinuity 
in the derivative dS/dT'—=o/T where o is the Thomson heat. Experi- 
mentally, however, it is difficult to gain any precise information by this 
means, for S is itself already the derivative of an experimentally measured. 


Fig. 6 
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The anomalous component of resistance as a function of temperature. 
A. ‘To test the validity of r4=f,(7¢) (6/T)”. 
B. To test the validity of r4=f,(r¢) exp (W4/RT). 
C. Showing the validity of r4—=f(r¢) (@-T)”. 
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e.m.f. Inspection of the values of the absolute thermoelectric power of 
alloys, Cu+-0-004 atomic % Fe and Cu+0-03 atomic °% Fe obtained from 
Borelius’ measurements (adjusted for the S values of the ‘ silver normal ’ 
_ against which measurements were made) shows that a discontinuous 
change in dS/d7' would not be impossible at the characteristic temperature 
expected from our measurements or a few degrees below it, but it does 
not appear possible to analyse the available data more critically. 

In conclusion let us reiterate that the anomalous scattering of very small 
quantities of iron in copper appears to be the most remarkable yet 
observed by us both as regards size and position of resistance minimum 
and variation of thermoelectric power. The size of the anomalous 
component at 4°K in the alloy containing 0-044 atomic °%% Fe is equal 
to the thermal scattering in pure copper at ~ 48°xK (fig. 4). That the 
effects are not due only to the presence of a transition metal which may 
have unfilled holes in the d shell is shown by Cu—Ni alloys in which no 
such striking effects are found, even though magnetic susceptibility 
measurements indicate that the ‘holes’ in the d shell of nickel are not 
filled in dilute alloys of nickel in copper. In addition, lead and bismuth 
in copper give anomalous effects which approach in magnitude those due 
to iron and with these solutes there can be no question of * holes’ in the 
d shell. We look forward with lively interest to the results of investi- 
gations of the magnetic susceptibility of these alloys at low temperatures 
which are being undertaken in our laboratory to throw more light on 
this problem. 
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Division of Physics, National Research Council, Ottawa* 
[Received March 12, 1955] 


ABSTRACT 


We have determined the temperature (7'min) at which the resistance 
minimum occurs for a large number of dilute copper alloys containing 
Fe, Ga, In, Si, Ge, Sn, Pb and Bi as solutes. Data on the absolute 
thermoelectric power and size of the resistance minimum of these alloys 
have previously been published. For alloys of Fe, Bi and Pb, we find 
that (T'min)" 0% rg where n ~ 5 and rg is the normal residual resistivity 
increase due to the added solute. The variation of 7min with solute 


concentration is much smaller for Ga, In, Ge and Sn as solutes than for 
Fe, Bi and Pb. 


$1. InrRopUCTION 


THE earliest investigations of Tin, the position at which the resistance 
minimum occurs, in relation to impurity content or residual resistivity 
(yes) were confined to gold (cf. de Haas and van den Berg 1936, 1937, 
Giauque, Stout and Clark 1937, Stout and Barieau 1939, MacDonald and 
Templeton 1951, MacDonald 1951). More recently the scope has been 
broadened to include in particular copper, silver and magnesium, as 
parent metals (cf. Gerritsen and Linde 1951, 1952, MacDonald 1952, 
Herlin 1953). These investigations have shown generally that T'min is 
a function of the specific impurity atom and its concentration. 

We have already examined the variation of the size of the resistance 
minimum with solute concentration in a series of dilute copper alloys 
with various solutes (MacDonald and Pearson 1955 a, b, Pearson 1955). 
In the present note, we give data on the temperature at which the 
resistance minimum occurs in these alloys. The methods of preparing 
and examining very dilute alloys of Fe, Ga, In, Si, Ge, Sn, Pb and Bi 
in copper have already been fully discussed. In the present work, 
when the minimum is of small size, it is often not possible to fix its 
position to better than +1-5°K. 
SE ee A ee 
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§2. EXPERIMENTAL RESULTS AND DISCUSSION 


The variation in position of the resistance minimum in Cu-Fe alloys 
as a function of iron concentration, shown in fig. 1, is the most striking 


Fig. 1 
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The temperature at which the resistance minimum occurs in dilute Cu—Fe 
alloys as a function of the residual resistance of the alloys. 


which we have found in a series of alloys. It changes from 13-5°K for 
0-002 atomic % Fe to ~ 28°K for 0-05 atomic °% Fe in solid solution and 


Fig. 2 
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‘The temperature at which the resistance minimum occurs as a function of ro, 
(the increase of residual resistance due to the normal ’ scattering caused 
by (1) bismuth, (2) iron in solid solution) showing rg & (7'min)”. 
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obeys rather exactly the law rg=C(T'min)”" where n ~ 5-2 (fig.2) and r¢ is 
the ‘normal’ impurity scattering due to the component which gives 
rise to the resistance minimum. The behaviour of lead and bismuth as 
solutes is similar to that of iron (e.g. see fig. 3). In the case of bismuth 
we find r¢=C(T min)” where n ~ 4:8 (fig. 2). 


Fig. 3 
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The temperature at which the resistance minimum occurs in dilute Cu-Bi 
alloys as a function of the residual resistance of the alloys. 


In alloys having Ga, In, Ge and Sn as solutes we find that 7'min increases 
rapidly from ~ 10-5°K for ~ 0-002 atomic % solute to a constant or very 
slowly rising value between 13° and 15°K in alloys containing up to 
~ 0-5 atomic °%, solute. The solute concentration at which 7Tmin becomes 
approximately constant coincides with the concentration at which the 
actual size of the minimum saturates. 

In the more dilute alloys containing Ga, In, Ge and Sn as solutes 
prepared under reducing conditions, the small amount of iron entering 
solid solution (~ 0-002 wt. °% present originally in the copper as an 
oxide) appears to control the position of the minimum which then lies 
generally between 13-5° and 15°K (over 90°, of the measurements lie in 
this range). 

The data shown above for alloys of iron and bismuth in copper now 
make it possible to examine critically a suggestion by MacDonald (1949, 
1951) that it might be possible to treat the resistance minimum formally 
in somewhat the same way as an applied magnetic field is introduced 
in the theory of magneto-resistance by writing 
Bae 
ro 
where rp is the resistive component due to ideal thermal scattering and. 
ry is the resistive component due to ‘ normal ’ impurity scattering. 


t=tet Tob 
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This led to the prediction that rmin & (Tmin)” where n has the value 
~ 5, given by the expression for thermal scattering of ideally pure 
copper (r=A.7™). We see that this prediction is rather well obeyed 
for bismuth and iron as solutes in copper. where rmin changes by a factor 
of 30. The model, however, does not lead quantitatively to the observed 
temperature dependence of the anomalous increase of resistivity of these 
alloys below Tin (cf. the previous paper) and it does not appear possible 
to attain agreement by any simple modification of the above expression. 
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CV. CORRESPONDENCE 


Gamma-rays from Inelastic Scattering of D-T’ Neutrons from 
Carbon 12 and Oxygen 16 


By L. E. Beautan, D. Hicks and B. MILMAN 
Clarendon Laboratory, Oxford 


[Received July 1, 1955] 


THE gamma-ray spectra resulting from the inelastic scattering of D-T 
neutrons by !2C and !6O have been investigated. Simple spectra con- 
sisting of strong ground state transitions from the 4-43 Mey level in #C 
(Thompson and Risser 1954, Battat and Graves 1955) and the 6-13 Mev 
level in 160 (Thompson and Risser 1954) have been reported previously. 
More complex spectra were however obtained by Scherrer e¢ al. (1953). 
The aim of the present investigation was to search for weaker radiations. 
from the higher excited states and in particular to establish the energy of 
the second excited state of 12C which has an astrophysical interest (Opic 
1952, Hoyle 1954). In fact, for both #C and 18O, only a single line was. 
observed and an upper limit is placed on the intensity of radiation from 
the second excited state of #C. 


EXPERIMENTAL ARRANGEMENT 


14 Mev neutrons were produced by bombarding a zirconium tritide 
target with 300 kev deuterons. Since gamma radiation is produced by the 
neutrons in the target assembly itself this was constructed as lightly as 
possible. It consisted of a thin aluminium can supported by a glass tube 
and cooling was provided by an air blast. 

A large scale ring geometry was employed in these experiments. The 
gamma radiation was observed in a direction perpendicular to the motion 
of the primary neutrons entering the scatterer thus eliminating the effects 
of Doppler shift (see fig. 1). 

The garnma radiation was detected with a 1 in. cube Nal(TI) crystal 
and a Dumont 6292 photomultiplier, and the resulting pulse height 
distribution was analysed with a 25 channel kicksorter. 


EXPERIMENTAL TECHNIQUE 


Pulse height distributions were recorded with and without the scatterer 
and were normalized to the same total neutron emission from the source. 
From these measurements the difference D was calculated. 

D (Counts per channel with scatterer) — 
— 1 . 
(Counts per channel without scatterer) 
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Fig. 1 
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(1) Glass tube. (2) Aluminium can. (3) Zirconium tritide target. 
(4) Scattering annulus. (5) Photomultiplier case. (6) Preamplifier. 


Fig. 2 
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A plot of the difference D will give the spectrum of the gamma-rays 
associated with inelastic scattering superimposed on a general background 
due to the scattered neutrons. 


MEASUREMENTS WITH CARBON 


A solid graphite annullus was used as scatterer in the C experiment. 
The spectrum shows the presence of only one gamma-ray of energy 
4-4 mev, leading to the typical one and two quantum escape pair peaks. 
This corresponds to the excitation of the first level in #C (see Ajzenberg 
and Lauritsen 1955). A search was made for a gamma-ray of 3-1 Mev 
which would result from a 7-65 mev level decaying by cascade through 
the 4:43 Mev state. 

However in this region of the spectrum the background is rising 
sharply, and a weak line might escape detection. An upper limit to the 
intensity of the 3-1 Mev line can be placed at 10°% of the intensity of the 
4-4 Mev gamma-ray. 

These results are consistent with the measurement of Battat and 
Graves (1955). However no evidence was found for the presence of 
gamma-rays of 2-8, 6-0 and 7-0 Mev reported by Scherrer, Theus and Faust 
(1953). 

MEASUREMENTS WITH OXYGEN 


Water, contained in a thin aluminium annullus was used as the 160 
scatterer. The spectrum obtained reveals the presence of a strong line of 
6-1 Mev energy. This is consistent with a transition from the second 
level in 1°O (6-13 Mev) to the ground state (see Ajzenberg and Lauritsen 
1955). No gamma transition is to be expected from the first level 
(6-05 Mev) to the ground state since this is a forbidden transition. However 
the resolution is not sufficient to exclude the possibility of some contri- 
bution to the observed peak by radiation from 6-9 and 7-1 Mev states 
which may also be excited. No evidence was found for gamma-rays of 
3-0, 3-8, and 5-2 Mev reported by Scherrer, Theus and Faust (1953) 
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